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ABSTRACT 


This  is  the  first  of  three  volumes.  Volumes  II  and  III  contain 
other  aspects  cf  the  study:  description  of  the  LITE  codes  and  the  two 
computer  programs  that  were  prepared  for  use  in  computing  cross  sections 
for  light  scattering  by  spherical  aerosol  particles  with  a  complex  in¬ 
dex  of  refraction. 

The  Monte  Carlo  codes  LITE-I  and  LITE-II  were  used  to  compute  light 
transmission  data  for  point  and  plane  sources.  The  results  of  these  cal¬ 
culations  showed  that  there  is  no  correlation  between  the  scattered 
transmission  data  for  point  isotropic  and  plane  parallel  sources. 

The  LITE  codes  were  used  to  analyze  experimental  data  on  light  trans¬ 
port  in  the  atmosphere.  Reasonably  good  agreement  was  obtained  for  those 
cases  where  data  were  available  to  adequately  describe  the  atmosphere  at 
the  time  of  the  experiments. 

Calculations  were  performed  to  determine  the  angular  dependence  of 
the  number  albedo  from  thick  cumulus  clouds  for  0.45p  light  incident  at 
various  angles  to  the  cloud.  The  polar  angular  distributions  of  the 
reflected  photons  were  found  to  be  cosine  distributions  and  the  total 
number  albedo  was  fitted  by  the  expression 

o(0  )  -  0.639  -  0.154cose 
0  o 

where  9q  is  the  angle  of  incidence. 

Studies  were  performed  to  determine  the  sensitivity  of  the  LITE-I 
calculated  scattered  intensities  for  a  point  isotropic  monochromatic 
source  to  changes  in  the  source  and  receiver  altitudes,  aerosol  number 
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density,  aerosol  phase  function,  altitude  variation  of  the  aerosol  scatter¬ 
ing  coefficient  and  the  altitude  of  the  bottom  of  a  cloud  layer  above  the 
source  point.  It  was  found  for  the  source-receiver  geometry  considered 
in  these  studies  that  the  LITE  calculations  were  more  sensitive  to  changes 
in  the  number  density  of  the  aerosol  particles  than  to  changes  in  the 
shape  of  the  aerosol  phase  function. 
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PREFACE 


During  the  period  1  August  1965  to  31  August  1966  Monte  Carlo  studies 
were  performed  to  determine  light  transport  in  the  atmosphere  under  vari¬ 
ous  environmental  conditions.  These  studies  consisted  of  1)  correlation 
analysis  of  light  transport  from  a  point  isotropic  source  and  a  plane 
parallel  source  to  determine  the  comparability  of  solar  light  transmission 
data  and  transmission  properties  for  thermal  radiation  nuclear  weapons, 

2)  development  of  machine  codes  for  calculation  of  phase  functions  and 
scattering  ano  absorption  coefficients  for  spherical-homogeneous  aerosol 
particles  with  a  complex  index  cf  refraction,  3)  an  analysis  of  experi¬ 
mental  field  a-Jta  on  light  transmission,  4)  parametric  studies  to  deter¬ 
mine  the  specific  influence  of  ground  and  cloud  albedo,  cloud  height  and 
aerosol  number  density  and  particle-size  distribution  on  the  transport  of 
light  in  the  atmosphere,  5)  modifications  to  the  LITE  codes  to  increase 
their  application  to  a  wider  range  of  atmospheric  transport  problems  and 
6)  the  development  of  a  machine  program  for  use  in  converting  the  scattered 
intensities  computed  by  the  LITE  codes  for  a  given  ground  albedo  to  data 
giving  scattered  intensities  and  scattered  fluxes  for  other  ground  albedos. 
The  results  of  these  studies  are  presented  in  this  report,  which  is  divided 
into  three  volumes.  The  first  volume  describes  the  results  of  items  1, 

3  and  4  outlined  above.  The  second  volume  describes  the  machine  programs 
developed  for  use  in  calculation  of  aerosol  cross  sections.  The  third 
volume  contains  utilization  instructions  for  the  modified  versions  of  the 
LITE  codes  and  for  the  code  developed  to  convert  the  LITE  results  to  data 
giving  scattered  intensities  and  fluxes  for  other  ground  albedos. 
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I .  INTRODUCTION 


The  present  report  is  Volume  I  of  a  three-volume  annual  report  pre¬ 
pared  under  Contract  No.  DA  28-043  AMC-00240(E) .  The  studies  described 
in  this  volume  were  performed  for  the  purpose  of  continuing  to  check  the 
accuracy  of  the  Monte  Carlo  codes,  LITE-I  and  LITE-II,  and  to  investigate 
the  importance  of  the  various  parameters  that  affect  the  transport  of 
light  in  the  atmosphere.  The  LITE  codes  that  were  used  in  these  studies 
are  described  in  Volume  III  of  this  report  and  in  Reference  1. 

The  first  study  undertaken  during  the  reporting  period  was  directed 
toward  determining  if  a  correlation  exists  between  atmospheric  trans¬ 
mission  data  for  point  and  plane  sources.  The  results  of  that  study  are 
given  in  Chapter  II. 

Although  there  is  currently  available  in  the  literature  a  large  body 
of  measured  data  on  light  transport  in  the  atmosphere,  most  of  that  data 
is  not  suitable  for  use  in  checking  out  the  results  of  theoretical  calcu¬ 
lations.  The  main  difficulty  with  the  available  measured  data  is  the 
almost  complete  lack  of  information  on  the  aerosol  number  density  and  size 
distribution  and  their  variation  with  altitude.  In  nany  cases  transmission 
measurements  were  taken  for  only  one  source  altitude  and  as  a  result,  it 
is  almost  impossible  to  infer  the  variation  of  the  extinction  cross  sec¬ 
tion  with  altitude  unless  one  can  describe  the  atmosphere  by  a  model,  such 
as  that  suggested  by  Elterman  (Reference  2). 

In  spite  of  the  many  uncertainties  in  the  atmospheric  conditions  for 
a  given  set  of  measurements',  several  sets  of  measured  data  were  selected 
for  analysis.  Chapter  III  presents  the  results  of  Monte  Carlo  calculations 
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performed  for  analysis  of  measured  data  from  several  experiments.  Section 
3.1  describes  an  analysis  of  transmission  data  reported  by  Cantor  and 
Petriw  (Reference  3)  for  point  isotropic  sources,  Section  3.2  describes 
an  analysis  of  some  of  the  measured  receiver  response  data  from  a  search¬ 
light  experiment  reported  by  Elterman  (Reference  4),  Section  3.3  discusses 
an  analysis  of  some  measured  sunlight  transmission  data  from  an  Air  Force 
Cambridge  Research  Laboratory  high  altitude  balloon  flight  and  Section  3.4 
shows  a  comparison  between  two  different  calculations  of  the  dependence 
of  light  transmitted  and  reflected  from  an  aerosol  atmosphere  on  the  shape 
of  the  aerosol  phase  function. 

The  need  for  information  on  the  angular  distribution  of  cloud  albedos 
for  use  in  the  LITE-1  code  resulted  in  a  study  to  determine  the  angular 
distribution  of  the  number  albedo  for  various  angles  of  incidence  of  0.45 
micron  light  to  a  1  km  thick  cumulus  cloud.  The  results  of  the  cloud 
albedo  studies  are  presented  in  Chapter  IV. 

Since  the  LITE-I  code  is  to  be  used  to  study  the  atmospheric  trans¬ 
port  of  the  thermal  radiation  emitted  by  a  nuclear  weapon,  it  is  of  con¬ 
siderable  importance  to  have  a  knowledge  of  how  variations  in  such  para¬ 
meters  as  the  source  altitude,  receiver  altitude,  cloud  height,  aerosol 
number  density,  aerosol  size  distribution  and  altitude  variation  of  the 
aerosol  number  density  affect  the  transmission  properties  of  the  atmosphere. 
Chapter  IV  describes  the  results  of  Monte  Carlo  calculations  that  were  per¬ 
formed  for  the  purpose  of  studying  the  influence  of  the  above  mentioned 
parameters  on  the  atmospheric  transmission  of  light  emitted  by  point 
isotropic  sources. 
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II.  CORRELATION  STUDIES  FOR  POINT  AND  PLANE  SOURCES 


The  availability  of  a  large  body  of  measured  data  (Reference  5)  on 
the  transmission  of  sunlight  through  the  atmosphere  leads  one  to  inquire 
if  this  data  could  be  used  to  compute  the  atmospheric  transmission  of  the 
thermal  radiation  produced  by  a  nuclear  detonation  in  the  atmosphere. 


The  scattered  light  intensity  at  a  receiver  depends  on  the  source- 
receiver  geometry,  the  variation  of  the  scattering  and  absorption  cross 
sections  with  altitude,  and  the  source  wave  length.  The  direct  radiation 
at  the  receiver  depends  only  on  the  number  of  mean-f ree-paths  between  the 
source  and  receiver.  The  direct  intensity  transmitted  through  an  atmos¬ 
phere  of  p  mean-free-paths  in  normal  thickness  for  a  broad  beam  monodirec- 
tional  source  of  one  photon  per  unit  area  incident  to  the  top  of  the 
atmosphere  at  angle  0o  is  given  by 


VV 


T  -psec0 
I  e  o 

o 


> 


-2 

where  the  incident  intensity  I  *  sec0Q  (photons  cm  ) .  For  a  point  iso¬ 
tropic  source  placed  at  the  top  of  the  atmosphere,  the  direct  intensity 
at  a  receiver  located  at  the  bottom  of  the  atmosphere  and  at  a  slant  dis¬ 
tance  R  from  the  source  is  given  by 


IT  (R  ■  hsecO  ) 
D  o 


-psec0 
e _ o 

2  2 
4ttIi  sec  0 


o 

where  h  is  the  normal  thickness  of  the  atmosphere  and  0q  is  the  angle 
between  a  normal  to  the  ground  through  the  source  and  the  source-receiver 


line.  It  is  seen  that 


I » 


I_(0  ) cos 6  -  4irh^sec^0  •  I^CR  ■  hsec0  )  .  (1) 

The  requirement  that  the  point  source  be  placed  at  the  top  of  the  atmos¬ 
phere  can  be  removed  when  calculating  the  direct  intensity  provided  that 
the  slant  mean-f ree-patn  distance  p'  between  the  source  and  receiver  be 

the  same  as  that  given  by  p'  =  psec0  where  0  is  the  angle  of  incidence 

o  o 

for  the  plane  source.  If  the  point  source  is  located  at  altitude  h^ 
ana  the  receiver  at  h^  and  the  angle  between  a  normal  to  the  ground 
through  the  receiver  and  the  source  receiver  line  is  a  then  the  mean- 
free-path  distance  between  the  source  and  receiver  is  given  by 

hr2 

p'  =  seen  ;  I(h)dh 
hl 

where  1(h)  is  the  extinction  cross  section  at  altitude  h.  The  slant  dis¬ 
tance  R  is  given  by 

R  =  (h^  -  h^)seca  .  (2) 

Therefore,  when  the  mean-f ree-path  distance  between  the  point  source  and 
receiver  is  the  same  as  that  along  a  slant  path  through  the  atmosphere 
for  a  plane  source,  we  can  estimate  the  direct  intensity  for  the  point 
source  by  use  of  the  relation 

COS0 

VR'p,) "  V6o’p,)  (3) 

4irR 

where  1^(0, p  * )  is  the  direct  intensity  at  a  depth  of  p'  mean-f ree-paths 
into  the  atmosphere  for  a  broad  beam  monodirectional  source  of  one  photon 
unit  area  incident  at  angle  to  the  top  of  the  atmosphere. 
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If  the  atmosphere  was  a  pure  absorbing  medium  so  that  no  scattering 
occurs,  then  measured  transmission  data  for  sunlight  incident  on  the  top 
of  the  atmosphere  could  be  used  to  predict  the  transmission  for  point 
sources.  The  mean-free-path  distance  from  altitude  h  to  the  top  of  tne 
atmosphere  could  be  obtained  from  measurements  of  sunlight  intensities  at 
various  altitudes. 

The  earth's  atmosphere  is  not  a  pure  absorbing  medium,  and  scattering 
is  an  important  process  by  which  light  is  transported  through  the  atmos¬ 
phere.  It  would  be  instructive  to  investigate  the  validity  of  Equation  1 
when  the  direct  intensities  are  replaced  by  the  scattered  intensities. 

In  addition  to  aerosol  and  Rayleigh  scattering,  another  important  scatter¬ 
ing  process  is  the  reflection  of  both  the  direct-beam  and  the  air-scattered 
light  by  the  ground  surface  and  clouds. 

The  scattered  light  intensity  transmitted  through  the  bottom  of  the 

-2 

atmosphere  for  a  plane  isotropic  source  emitting  0.5  photons  cm  into 
the  atmosphere  may  be  computed  from  a  knowledge  of  Ig(9o),  the  scattered 
intensity  transmitted  through  the  bottom  of  th«^  atmosphere  for  a  plane 
parallel  source,  and  I  (R  =  hsec0  ),  the  scattered  intensity  transmitted 

b  O 

through  the  bottom  of  the  atmosphere  at  a  slant  distance  R  from  a  point 
isotropic  source  located  at  the  top  of  ths  atmosphere. 

The  transmitted  intensity  for  a  plane  isotropic  source  is  given  by 

w/2  I-(8  )  tt/2 

I.  *  2  it  /  sine  de  -  ”  /  Ic(0  )sine  de  , 

S  >  4tt  o  o  l  ‘  S'  o  O  o 


the  equation 
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where  7-  is  the  source  strength  per  steredian  in  direction  6  and  I _ ( e  ) 

4TT  0  SO 

is  the  scattered  intensity  transmitted  through  the  bottom  of  the  atmosphere 
-2 

per  photon  cm  incident  at  direction  @o  to  the  top  of  the  atmosphere. 

The  scattered  intensity  transmitted  through  the  bottom  of  the  atmos- 

-2 

phere  for  a  plane  isotropic  source  emitting  0.5  photons  cm  into  the 
atmosphere  is  also  given  by  Integrating  the  transmitted  scattered  intensity 
for  a  point  isotropic  source  located  at  the  top  of  the  atmosphere  over  the 
area  of  the  bottom  of  the  atmosphere.  If  h  is  the  thickness  of  the  atmos¬ 
phere,  R  the  slant  range  between  the  source  and  receiver,  and  6  is  the 

o 

angle  between  R  and  a  normal  to  the  atmosphere  through  the  source,  then 

*/2  2  12 
T  *  2tt  /  i_(R  ■  hsec0  )R  tan0  d6  =7  /4ttR  I_(R  =  hsec0  )tan0  d0 
S  0S  0  oo2''  S  000 

To  show  that  the  above  integrals  are  equal,  LITE-I  calculations  were  run 

to  determine  Ig(R)  for  a  point  isotropic  source  located  at  the  top  of  a 

0.5  mean-free-path  thick  Rayleigh  atmosphere  and  for  receivers  placed 

parallel  to  the  ground  at  ground  level.  Data  for  I„(0  )  were  obtained 

o  o 

from  integration  of  the  differential  scattered  intensity  data  reported 
by  Coulson  et_  al.  (Reference  6).  The  results  of  the  integrations  are 
shown  below  for  cases  where  the  groi nd  albedo  was  0.0  and  0.8. 

Scattered  Intensity  Transmitted  Through  a  0.5  Mean- 
free-path  Thick  Rayleigh  Atmosphere  for  a  Plane  Isotropic  Source 

-2  -2 
(photons  cm  /0.5  source  photons  cm  ) 


Transmitted 

Intensity 

Ground 

From  Coulson 

From  LITE-I 

Albedo 

Data 

Calculations 

0.0 

0.332 

0.351 

0.8 

0.561 

0.568 
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It  is  seen  that  integrations  over  the  plane  and  point  source  data  give 
identical  results  within  the  accuracy  cf  the  graphical  integration  method 
used.  The  fact  that  the  two  integrals  are  equai,  that  is, 

,  tt/2  ,  tt/2  » 

4  /  Ic(e  )sinO  d 6  -  4  /  4nR  (e  ) I (R  =  hsece  )tan0  d6  , 

2gbO  O  O  t  0  o  oo 

and  the  fact  that  the  direct  intensities  from  the  two  sources  are  corre¬ 
lated,  leads  one  to  inquire  if  the  scattered  intensities  from  the  two 
sources  are  also  correlated. 

If  the  two  integrals  were  equal  for  all  values  of  the  upper  limit 
of  the  interval  of  integration  less  than  tt/2,  then  differentiation  of 
both  integrals  would  give 

I0 (0  )sin0  *  4ttR2(0  )Ic(R  -  hsec6  )tar.0  , 

Soo  oS  oo 

and  a  correlation  would  be  obtained  for  point  and  plane  source  data.  The 
fact  that  these  integrals  are  equal  only  for  the  special  case  where  the 
upper  limit  of  the  integration  interval  is  tt/2  is  shown  in  Figure  1.  It 
is  seen  from  the  two  curves  shown  in  Figure  1  that 

e  0 

/  Ic(6  )d(cos6  )  if  f  4ttR  (0  )sec0  IC(R  *  hsec0  )d(cos0  ) 

0  b  O  O  q  0  O  D  o  o 

2 

for  values  of  0<tt/2.  It  is  also  seen  in  Figure  1  that  4ttR  sec©  I(R=hsec0  ) 

o  o 

and  Ic(0  )  are  equal  only  at  two  values  of  cos©  .  Therefore,  it  can  be 
o  o  o 

stated  that  no  real  correlation  exists  between  the  scattered  intensities 
for  plane  parallel  and  point  isotropic  sources. 

A  series  of  LITE-I  and  LITE-II  problems  were  run  for  0.45  and  0.65 
micron  plane  parallel  and  point  isotropic  sources  in  an  Elterman  model 
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atmosphere  (Reference  2).  The  scattered  transmission  data  for  each  source 

were  plotted  as  a  function  of  the  angle  0q  and  the  results  obtained  showed 

that  there  was  no  correlation  between  1.(9  )  and  IC(R  -  hsece  )  for  both 

so  b  o 

source  wave  lengths  studied. 

As  a  result  of  the  analysis  described  above,  it  must  be  concluded  that 
significant  errors  would  be  introduced  by  using  measured  sunlight  trans¬ 
mission  data  to  approximate  the  atmospheric  transmission  of  the  thermal 
radiation  produced  by  nuclear  weapon  detonations  in  the  atmosphere. 
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III.  COMPARISON  OF  CALCULATED  AND  MEASURED  TRANSMISSION  DATA 

A  measure  of  the  accuracy  of  any  machine  procedure  for  calculating 
the  scattering  of  light  by  the  atmosphere  is  its  ability  to  predict  mea¬ 
surements.  As  a  part  of  the  continuing  effort  to  validate  the  LITE  codes, 
several  sets  of  experimental  data  were  selected  for  analysis.  Discussions 
of  the  LITE  calculations  and  their  comparison  with  experimental  data  are 
given  in  the  following  sections. 

One  of  the  difficulties  that  arise  when  performing  an  analysis  of 
measured  data  is  the  almost  complete  lack  of  information  concerning  the 
atmospheric  parameters  that  describe  the  atmosphere  at  the  time  the  mea¬ 
surements  were  taken.  In  an  attempt  to  overcome  this  difficulty,  com¬ 
parisons  were  also  made  with  the  results  of  calculations  produced  by 
other  machine  codes. 

3. 1  Analysis  of  Measured  Data  for  Point  Sources 

Calculations  were  run  using  the  LITE-I  code  to  predict  the  scatter¬ 
ing  of  yellow  light  (A  =  0.55  microns)  in  air  over  a  snow-covered  terrain 
for  comparison  with  measurements  (Reference  3)  taken  at  Camp  Haven, 
Wisconsin,  on  30  and  31  March  19o5.  The  source  was  located  in  the  tail 
of  an  airplane  flying  at  an  altitude  of  2650  ft  and  the  receiver  was 
located  at  an  altitude  of  20  ft. 

The  measured  direct  beam  was  used  to  obtain  an  estimate  of  the  ex¬ 
tinction  coefficient  for  each  day.  The  aerosol  and  extinction  coeffi¬ 
cients  as  given  by  Elterman's  model  were  modified  in  order  to  duplicate 
the  measured  extinction  coefficients  for  each  day. 
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The  measured  data  for  30  March  gave  a  value  of  0.15  mile  *  for  the 
extinction  coefficient.  The  measured  extinction  coefficient  for  31  March 
was  0.11  mile  \  The  atmosphere  was  clear  on  the  night  of  30  March  and 
the  surface  albedo  was  estimated  to  be  about  50  percent.  There  was  an 
overcast  at  3000  ft  altitude  on  the  night  of  31  March  and  the  surface 
albedo  was  estimated  to  be  about  25°.  The  phase  function  for  light 
scattering  at  ground  level  was  measured  on  both  evenings.  The  calculations 
for  31  March  assumed  that  the  overcast  was  a  cloud  of  at  least  1  km  in 
thickness  with  the  bottom  of  the  cloud  at  a  height  of  3000  ft.  Calculated 
fluxes  for  both  evenings  were  obtained  for  a  flat  plate  receiver  which 
was  assumed  to  be  normal  to  the  source-receiver  axis. 

The  calculated  and  measured  flux  data  for  30  March  are  compared  in 
Figure  2.  The  measured  data  giving  the  total  transmittance  (receiver 
response  times  slant  distance  squared)  were  arbitrarily  normalized  to 
give  a  value  of  1.0  at  a  slant  range  of  R  =  2.25  miles.  Therefore,  the 
calculated  total  transmittance  data  for  a  ground  albedo  of  0.5  were 
normalized  to  the  measured  data  at  R  *  2.25  miles.  The  estimated  albedo 
for  the  ground  on  the  day  the  measurements  were  taken  was  0.5.  Although 
there  are  large  fluctuations  in  the  measured  data,  it  is  seen  that  the 
falloff  of  the  total  transmittance  with  slant  distance  as  given  by  the 
calculated  data  is  in  good  agreement  with  the  measured  data.  The  large 
fluctuations  in  the  measured  data  make  it  impossible  to  estimate  the 
value  of  the  ground  albedo  from  a  comparison  of  the  measured  and  cal¬ 


culated  data. 


NORMALIZED  MEASURED  RESPONSE  TIMES  S  \NT  RANGE  SQUARED 
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Mg.  ?. .  Comparison  of  Calculated  and  Measured  Receiver  Response  for  AfmospheMc  Transmission 
of  Yellow  Light  in  a  Wisconsin  Area:  30  March  1^65 


Figure  3  presents  a  comparison  of  the  calculated  and  measured  total 
transmittance  data  for  yellow-light  scattering  on  March  31,  1965.  The 
calculated  and  measured  data  giving  the  total  transmittance  were  normalized 
to  a  value  of  1  at  R  ■  2.25  miles.  It  is  noted  that  there  is  a  wider 
spread  between  the  calculated  total  and  direct  transmittance  curves  in 
Figure  3  for  March  31  than  that  shown  in  Figure  2  for  March  30.  This 
difference  results  from  the  presence  of  the  cloud  located  at  a  neight 
of  3000  ft  for  the  case  shown  in  Figure  3.  The  comparison  between  the 
measured  and  calculated  values  of  the  total  transmittance  vs  R  indicates 
that  the  calculated  data  for  a  ground  albedo  of  0.25  are  about  15  per¬ 
cent  low  for  slant  ranges  between  4  and  15  miles. 

3.2  Analysis  of  Searchlight  Scattering  Data 

Elterman  (Reference  4)  has  carried  out  scattering  measurements  from 
a  searchlight  beam  to  determine  the  aerosol  properties  of  the  atmosphere. 

He  was  able  to  derive  the  variation  of  the  aerosol  scattering  coeffi¬ 
cient  with  altitude  from  the  measured  receiver  response  data  taking  into 
consideration  the  contribution  of  single  order  aerosol  and  molecular 
scattering  to  the  measured  response.  The  measured  response  data  obtained 
by  Elterman  for  8  May  1964  and  the  aerosol  cross  section  variation  with 
altitude  that  he  obtained  from  that  data  have  been  analyzed  by  use  of 
LITE- I. 

The  searchlight  scene  geometry  as  used  in  the  LITE-I  calculations 
is  shown  in  Figure  4.  The  searchlight  beam  in  the  light  scattering  cal¬ 
culations  was  defined  in  terms  of  a  9,  <f>  coordinate  system  so  that  it 
subtended  the  same  solid  angle  as  that  used  in  the  searchlight  experiment. 


NORMALIZED  MEASURED  RESPONSE  TIMES  SLANT  RANGE  SQUARED 
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R,  SLANT  RANGE  (miles) 


Fig.  3.  Comparison  of  Colculated  and  Measured  Receiver  Response  for  Atmospheric  Transmission 
af  Yellow  Light  in  a  Wisconsin  Areo:  31  March  1965 
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The  scattered  Intensities  at  the  receiver  were  recorded  for  those  pho¬ 
tons  reaching  the  receiver  through  solid  angles  defined  by  46  intervals 
in  the  receiver  polar  angle  about  the  source-receiver  line  and  6.8° 
intervals  of  the  azimuthal  angle  $  (see  Figure  5).  Although  At  was  not 
possible  with  the  ip,  b  coordinate  system  to  define  a  receiver  solid 
angle  that  did  not  vary  with  the  receiver  angle  of  elevation,  the  solid 
angle  formed  by  the  ret  iiver  polar  and  azimuthal  angular  intervals  used 
in  the  calculations  always  included  the  searchlight  beam.  Therefore, 
the  receivci  solid  angle  for  large  receiver  angles  of  elevation  actually 
viewed  a  scattering  volume  that  was  slightly  larger  than  the  volume 
that  would  have  been  formed  by  the  intersection  of  the  searchlight  beam 
and  a  constant  receiver  solid  angle  similar  to  that  used  in  the  experi¬ 
ment. 

Figure  6  shows  the  altitude  variation  of  the  various  cross  sections 
for  0.55  micron  wave  length  light  as  used  in  the  LITE-I  calculations. 

The  extinction  coefficient  was  taken  to  be  the  sum  of  the  aerosol  scat¬ 
tering,  Rayleigh  scattering,  and  ozone  absorption  coefficients.  The 
Rayleigh  scattering  and  ozone  absorption  coefficients  shown  in  Figure  6 
were  taken  from  a  tabulation  reported  by  Elterman  (Reference  2)  for  a 
"clear  standard  atmosphere."  One  of  the  aerosol  scattering  coefficient 
curves  shown  in  Figure  6  were  taken  from  Reference  2  and  the  other  was 
that  obtained  by  Elterman  from  an  analysis  of  the  searchlight  data  for 
8  May  1964. 

In  his  analysis  of  the  measured  data  Elterman  used  an  aerosol  phase 
function  that  had  been  obtained  from  measurements  reported  by  Reeger 
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CROSS  SECTION  fkm- 


18 


Fig.  6.  Altitude  Variation  of  the  Rcyleigh  and  Aerosol  Scottering  Coefficients  and  the  Ozone  Absorption 
Coefficient 


19 


and  Siedentopf  (Reference  7).  A  portion  of  the  normalized  Reeger  and 
Siedentopf  aerosol  phase  function  for  scattering  angles  between  75  and 
145  degrees  is  presented  in  Figure  7.  Also  shown  is  the  normalized 
phase  function  for  Rayleigh  scattering  in  the  same  scattering  angle 
range.  The  phase  functions  shown  in  Figure  7  have  been  normalized  so 
that  the  integral  of  the  phase  function  over  solid  angle  is  unity. 

Calculations  were  run  on  LITE-I  to  determine  the  single  scattered 
intensity  as  a  function  of  the  receiver  angle  of  elevation  that  resulted 
from  Rayleigh  scattering  only  and  from  both  Rayleigh  and  aerosol  scat¬ 
tering  where  the  aerosol  scattering  coefficient  variation  with  altitude 
was  that  obtained  by  Elterman  from  the  measured  searchlight  data.  The 
atmosphere  used  in  the  LITE-1  calculations  did  not  contain  any  ozone 
since  Elterman' s  earlier  analysis  was  based  only  on  a  Rayleigh  plus 
aerosol  scattering  atmosphere.  The  results  of  the  Monte  Carlo  calcula¬ 
tions  are  presented  in  Figure  8  where  they  can  be  compared  with  the 
shape  of  the  measured  receiver  response  data  as  a  function  of  the  re¬ 
ceiver  angle  of  elevation.  There  is  reasonably  good  agreement  between 
the  measured  and  calculated  data  when  the  calculated  data  was  normalized 
to  the  measured  data  at  0  *  54°.  The  largest  differences  occur  at  re¬ 
ceiver  angles  of  elevation  less  than  14  degrees. 

To  study  the  effect  of  neglecting  the  attenuation  of  the  scattered 
radiation  by  ozone  absorption  in  Elterman 's  calculation  of  the  aerosol 
coefficient,  another  problem  was  run  in  which  the  extinction  coefficient 
was  assumed  to  be  made  up  of  the  sum  of  the  Rayleigh  and  aerosol  scatter¬ 
ing  and  ozone  absorption  coefficients.  The  aerosol  scattering  coefficients 
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Fig.  7.  Phase  Functions  for  Rayleigh  Scattering  ond  Aerosol  Scotterina 


RECEIVER  RESPONSE 
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t),  RECEIVER  ELEVATION  ANGLE  (deg) 


Fig.  8.  Comparison  af  Single  Scattering  Calculations  with  Measured  Data:  Na  Ozone  Absorption 
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used  in  the  calculations  were  those  obtained  from  the  measured  search¬ 
light  response  data. 

The  calculated  results  for  this  case  are  presented  in  Figure  9 
where  they  are  compared  with  the  shape  of  the  measured  receiver  response. 
The  calculated  single  scattered  intensity  summed  over  Rayleigh  and  aero¬ 
sol  scattering  were  normalized  in  Figure  9  to  the  measured  response  data 
for  a  receiver  angle  of  elevation  of  54°.  It  is  seen  that  the  total 
single  scattered  intensity  (Rayleigh  plus  aerosol)  falls  off  a  slight 
bit  faster  with  the  angle  of  elevation  than  does  the  measured  response 
data,  indicating  that  there  are  probably  some  multiple  scattering  ef¬ 
fects.  The  differences  between  the  total  (Rayleigh  plus  aerosol)  single 
scattered  intensity  and  the  single  Rayleigh  scattered  intensity  become 
smaller  with  increasing  angle  of  elevation,  which  indicate  that  the 
total  single-scattered  intensity  results  mainly  from  Rayleigh  scattering 
at  angles  of  elevation  greater  than  about  40  degrees. 

The  LITE-I  problem  described  above  was  rerun  to  obtain  the  total 
scattered  intensity  at  the  receiver  for  orders  of  scattering  up  to  20. 
Figure  10  shows,  as  a  function  of  the  receiver  angle  of  elevation,  a 
comparison  of  the  computed  total  scattered  intensity  with  the  shape  cf 
the  measured  response  data,  where  the  total  calculated  data  have  been 
normalized  to  the  measured  response  data  at  8  =  54  degrees.  It  is  seen 
that  the  multiple  scattered  intensity  variation  with  receiver  angle  of 
elevation  reproduces  very  well  the  shape  of  the  measured  response  data. 

It  is  also  seen  from  comparison  of  the  total  single  scattered  intensity 
curve  in  Figure  9  with  the  multiple  scattered  intensity  curve  in  Figure  10 


RECEIVER  RESPONSE 
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8,  RECEIVER  ELEVATION  ANGLE  (deg) 


Fig.  9,  Comparison  of  Single  Scaffering  Calculations  with  Measured  Dcta:  Orone  Absorption  Included 


RECEIVER  RESPONSE 


0,  RECEIVER  ELEVATION  ANGLE  (deg) 


Fig.  10.  Comparison  of  Multiple  Scattering  Calculations  with  Measured  Date:  Orane  Absorption 
Included 
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that  the  inctease  in  the  scattered  intensity  for  a  given  receiver  angle 
resulting  from  orders  of  scattering  greater  than  one  is  small.  The 
importance  of  orders  of  scattering  greater  than  one  increases  with  the 
receiver  angle  of  elevation.  At  an  angle  of  elevation  of  54  degrees 
orders  of  scattering  greater  than  one  contribute  about  12  percent  to 
the  calculated  total  intensity  whereas  at  a  receiver  angle  of  2°,  orders 
of  scattering  greater  than  one  contribute  only  about  two  percent  of  the 
total  scattered  intensity. 

From  a  comparison  of  the  calculated  and  measured  response  data 
shown  in  both  Figures  9  and  10,  it  appears  that  the  ozone  absorption 
and  multiple  scattering  effects,  which  were  not  considered  by  Elterman 
in  his  analysis  of  the  measured  data,  cancel  out  each  other.  The  good 
agreement  between  the  calculated  and  measured  date  as  seen  in  Figure  10 
indicates  that  the  aeroso.l  scattering  coefficient  variation  with  altitude 
is  accurately  given  by  Elterman' s  analysis  of  the  measured  response  data 
from  the  searchlight  experiment. 

The  altitude  variation  of  the  aerosol  scattering  coefficient  ob¬ 
tained  by  Elterman  from  the  searchlight  experiment  ?'s,  as  shown  in 
Figure  6,  considerably  different  from  that  given  by  his  "clear  standard 
atmosphere"  model.  To  investigate  the  effect  that  the  altitude  variation 
of  the  aerosol  scattering  coefficient  has  on  the  calculated  receiver 
response,  a  LITE-I  problem  was  run  to  determine  the  variation  of  the 
single  scattered  intensities  with  receiver  angle  of  elevation  for  an 
atmosphere  described  by  Elterman' s  model  for  a  "clear  standard  atmos¬ 
phere".  The  Reeger  and  Siedentopf  phase  function  was  used  to  describe 
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the  angular  distribution  for  aerosol  scattering.  The  results  of  the 
LITE-I  single  scattering  problem  for  the  "clear  standard  atmosphere" 
are  compared  in  Figure  11  with  the  results  obtained  for  an  atmosphere 
where  the  aerosol  scattering  coefficient  was  taken  to  be  that  given  by 
an  analysis  of  the  measured  response  data  from  the  searchlight  experiment. 
The  single  scattered  intensities  from  the  calculations  for  the  two 
different  aerosol  atmospheres  do  not  differ  significantly  except  in  the 
region  between  receiver  angles  of  elevation  of  two  to  ten  degrees.  For 
receiver  angles  of  elevation  above  ten  degrees  the  "clear  standard 
atmosphere"  has  a  smaller  aerosol  cross  section  than  does  the  atmosphere 
at  the  time  of  the  searchlight  experiment.  This  results  in  higher  single 
Rayleigh  scattered  intensities  due  to  the  decreased  attenuation  along 
the  path  from  the  source  to-scattering  volume-to-receiver . 

The  results  of  the  LITE-I  calculations  for  analysis  of  the  search¬ 
light  experiment  indicate  that  the  LITE-I  code  is  capable  of  predicting 
with  reasonable  accuracy  the  measured  receiver  response  data. 

3.3  Analysis  of  Measured  Sunlight  Transmission  Data 

The  LITE-II  code  was  used  to  analyze  the  results  of  a  high  altitude 
balloon  flight  carried  out  by  the  Air  Force  Cambridge  Research  Laboratory 
(Reference  8)  which  measured  various  brightness  parameters  at  100,000 
altitude.  The  experimental  data  gave  for  the  100,000  foot  level  the 
illuminance  (i.e.,  the  incident  radiation  for  a  photopic  spectral  sen¬ 
sitivity)  on  a  vertical  surface  as  a  function  of  the  solar  zenith  angle. 
Data  were  taken  for  cases  where  the  normal  to  the  vertical  surface  was 
facing  the  sun  (4>  *  0°)  or  facing  opposite  to  the  sun  direction. 


SCATTERED  INTENSITY  IN  RECEIVER  SOLID  ANGLE  (photons  cm'Vsource  photon) 
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Fig,  II.  Comparison  of  Single  Scattering  Calculations  for  Two  Different  Altitude  Distributions  of  the 
Aerosol  Coefficient 
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The  calculations  were  made  for  0.5  micron  wave  length  plane  parallel 
sources  since  that  wave  length  corresponds  approximately  to  the  peak  of 
the  photopic  curve.  The  atmosphere  was  described  by  the  Elterman  model 
for  a  "clear  standard  atmosphere".  The  phase  function  for  aerosol  scat¬ 
tering  was  described  by  Fraser's  (Reference  9)  Model  C  aerosol  size  dis¬ 
tribution.  The  ground  surface  reflectivity  was  estimated  to  be  about 
ten  percent,  but  during  the  balloon  flight,  there  was  a  change  in  cloud 
cover  which  could  have  changed  the  reflectivity  to  about  50  percent. 

The  LITE-II  output  gave  the  polar  and  azimuthal  angular  distribu¬ 
tions  of  the  scattered  intensities  at  an  altitude  of  30  km.  By  use  of 
the  Albedo  Conversion  Code  (Reference  10),  the  LITE-II  output  was  con¬ 
verted  from  intensities  to  fluxes  across  vertical  flat  plate  receivers 
oriented  so  that  <  -y  were  either  facing  the  sun  or  facing  opposite  to 
the  sun  direction.  The  scattered  flux  data  for  each  receiver  azimuthal 
direction  and  sun  zenith  angle  were  multiplied  by  the  secant  of  the 

solar  zenith  angle  to  give  fluxes  for  a  source  strength  of  one  photon 
-2 

cm  passing  through  a  surface  normal  to  the  sun  direction.  The  re¬ 
sults  so  obtained  for  ground  albedos  of  0.0,  0,5  and  1.0  are  presented 

in  Table  I.  The  direct  flux  incident  on  the  vertical  receiver  per 
2 

photon  per  cm  of  area  normal  to  the  zenith  direction  is  given  by  the 
equation 

N  (0  )  =  tan6  e  psec®o 
Do  1 

where  p  is  the  normal  mean-free-path  distance  from  an  altitude  of  30  km 
to  the  top  of  the  atmosphere.  For  0.5  micron  wave  length  light  p  =  0.004. 
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Table 

I.  Number  Current  Passing  Through 
at  an  Altitude  of  30  km:  A  ■  0. 

a  Vertical 

5  microns 

Surface 

(photon3  cm 

■2 

/source  photon  cm 

^  for  source  area  normal  to 

zenith  direction) 

Zenith 

Angle 

Direct 

Flux 

Scattered  Flux 

0Q(deg) 

4  -  0? 

Albedo 

4  =  0° 

4  -  180° 

30 

0.5774 

0.0 

9.13-2 

1.01-2 

0.5 

2.75-1 

2.96-1 

1.0 

4.99-1 

5.27-1 

40 

0.8391 

0.0 

1.30-1 

1.25-1 

0.5 

3.34-1 

3.42-1 

1.0 

5.83-1 

6.05-1 

50 

1.192 

0.0 

1.98-1 

1.92-1 

0.5 

4.38-1 

4.29-1 

1.0 

7.28-1 

7.19-1 

60 

1.732 

0.0 

3.31-1 

2.78-1 

0.5 

6.26-1 

5.79-1 

1.0 

9.85-1 

9.44-1 

30 
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Therefore,  very  little  error  would  be  introduced  by  taking  N^(6o)  *  tan6Q. 
From  the  data  given  in  Table  I  it  is  seen  that  although  the  scattered  flux 
incident  on  the  receiver  is  highly  dependent  on  the  ground  albedo,  the 
scfttered  flux  passing  through  either  side  of  the  flat  plate  receiver  is 
approximately  a  constant  for  a  given  source  zenith  angle  and  ground  albedo. 

A  comparison  of  the  calculated  data  for  a  ground  albedo  of  0.0  with 
the  measurements  is  shown  in  Figure  12.  The  calculated  total  fluxes  have 
been  normalized  to  the  measured  data  for  $  =  0°  at  a  zenith  angle  of  30° 
The  calculated  data  indicates  that  the  effective  ground  albedo  for  the 
measurement  was  close  to  zero.  The  measured  flux  for  <t>  =  180°  was  over¬ 
predicted  by  a  factor  of  1.4  at  6  =  30°  and  a  factor  of  2  at  6  =  60°. 

o  o 

The  overprediction  of  the  scattered  flux  for  $  =  180°  is  believed  to  be 
due  to  the  fact  that  although  the  peak  in  the  solar  spectrum  is  at  about 
0.5  microns,  nearly  76.5  percent  of  the  incident  energy  has  wave  lengths 
above  0.5  microns.  Since  the  scattering  cross  section  decreases  with 
increasing  wave  length,  it  is  expected  that  the  scattered  flux  would 
also  decrease  with  increasing  wave  length  provided  the  ground  albedo  was 
negligible.  The  wave  length  which  divides  the  solar  spectrum  into  two 
intervals  containing  one  half  of  the  incident  energy  is  0.72  microns. 

A  review  of  Couison,  Dave  and  Sekera's  (Reference  6)  calculations  of  the 
number  of  photons  reflected  from  a  Rayleigh  atmosphere  showed  that  for 
a  ground  albedo  of  0.0,  the  number  reflected  is  dependent  on  the  mean- 
free-path  thickness  of  the  atmosphere  and  the  angle  of  incidence.  An 
assumption  that  the  ext. action  cross  section  is  made  up  entirely  of 
Rayleigh  scattering  for  wave  lengths  of  0.5  and  0.72  microns  would 
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Fig.  12.  Comporison  of  Colculoted  ond  Meosured  llluminotion  os  seen  by  o  Verticol  Surface  ot 
on  Altitude  of  30  Km 
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result  in  about  1.5  times  as  many  photons  being  reflected  with  a  wave  length 
of  0.5  microns  as  will  be  reflected  with  a  wave  length  of  0.72  microns.  A 
significant  portion  of  the  solar  energy  incident  on  top  of  the  atmosphere 
with  wave  lengths  abcve  0.7  microns  will  be  absorbed  by  the  water  vapor 
and  carbon  dioxide  in  the  atmosphere.  Therefore,  it  is  reasonable  to 
expect  that  the  scattered  fluxes  computed  for  0.5  micron  wave  length 
light  would  be  higher  than  those  measured  during  the  balloon  experiment. 

The  percent  transmission  of  light  through  the  bottom  of  the  atmosphere 

as  a  function  of  the  source  angle  of  elevation  was  obtained  from  the  Monte 

Carlo  calculations  described  above  for  0.5  micron  wave  length  light  and 

from  Monte  Carlo  calculations  for  0.45  and  0.65  micron  wave  length  light 

The  atmosphere  for  the  0.45  and  0.65  micron  wave  length  sources  was 

assumed  to  be  that  given  by  Elterman  and  the  aerosol  phase  function  was 

-4 

taken  from  Mie  calculations  for  a  size  distribution  proportional  to  r 
The  results  for  a  ground  albedo  of  C.C  are  plotted  in  Figure  13  as  a 
function  of  the  source  elevation  angle.  It  is  seen  that  the  percent 
transmission  increases  with  the  source  wave  length  and  with  the  source 
elevation  angle.  The  average  percent  transmission  of  sunlight  for  vari¬ 
ous  solar  angles  of  elevation  was  obtained  from  an  analysis  of  measure¬ 
ments  reported  by  18  different  weather  stations  (Reference  5)  The  re¬ 
sults  are  compared  in  Figure  13  with  the  calculated  data  for  0.45,  0,50 
and  0.65u  wave  length  light.  The  measured  sunlight  transmission  data 
were  taken  from  data  compiled  in  Reference  5  for  0  to  30  percent  cloud 
cover.  There  was  no  information  available  on  the  value  of  the  ground 
albedo  for  each  reporting  station.  It  is  difficult  to  draw  any  definite 
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Fig,  i  3.  Transmissivity  as  a  Function  of  Solar  Elevation  Angle:  Ground  Aibedo  0,0 
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conclusions  about  the  agreement  between  the  measured  and  calculated  trans¬ 
mission  data  considering  that  the  atmospheric  conditions  varied  daily  at 
each  of  the  weather  stations.  It  should  be  pointed  out  that  the  trans¬ 
mission  for  wave  lengths  above  0.7  microns  should  be  less  than  that  com¬ 
puted  for  a  wave  length  of  0.65  microns  since  absorption  by  water  vapor 
and  C02  becomes  important  for  wave  lengths  above  0.7  microns. 

3.4  Comparison  of  LITE-II  Calculations  with  Irvine's  Calculations 

Irvine  (Reference  11)  has  conducted  a  study  to  determine  the  effect 
of  varying  the  shape  of  the  Mie  scattering  phase  function  on  the  angular 
intensities  transmitted  and  reflected  from  both  0.5  and  1  mean-free- 
path  length  th^ck  slab  atmospheres.  The  phase  functions  used  in  Irvine's 
calculations  sere  described  by  the  expression 


$ ( cose.)  «=  b$^G(g^, cot?a)  +  (l-b)<j>HG(g2,cosa) 


where  $K(,(g,cosa) 


(1+g  -2gcosa) 


3/2  • 


We  denote  as  case  D  the  results  obtained  using  valuesof  0.9724,0.824  and 
-0.55  for  b^,  g^  and  g^,  respectively  in  the  equation  for  4>(cosc:) .  The 
results  obtained  when  b^=  0.95,  g^a  0.9  and  g2=  0.75  are  denoted  as 
case  H.  The  phase  function  for  case  H  is  peaked  more  highly  forward 
than  is  the  phase  function  for  case  E.  Figure  14  shows  the  cumulative 
distribution  functions  for  cases  E  and  H. 


Two  LITE-II  problems  were  run  for  a  plane  source  of  light  incident 
normal  to  an  atmosphere  one  mean-f ree-path  in  thickness.  In  the  first 
problem  the  phase  function  for  case  E  was  used,  and  in  the  second 
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problem,  the  phase  function  for  case  H  was  used.  The  results  of  the  cal¬ 
culations  for  case  E  are  compared  in  Figure  15.  The  LITE-II  results  and 
Irvine's  results  for  case  H  are  compared  in  Figure  16.  The  transmit  ed 
angular  intensities  obtained  from  the  LITE-II  results  are  in  good  agree¬ 
ment  with  Irvine's  calculations.  The  differences  between  the  results 
obtained  from  the  LITE-II  calculations  and  Irvine's  calculations  for  the 
reflected  intensities  are  probably  due  to  poor  statistics,  since  less 
than  seven  percent  of  the  particles  undergoing  first  collisions  in  the 
atmosphere  scattered  through  an  agnle  greater  than  90®.  It  is  felt  that 
better  agreement  in  the  reflected  distribution  would  have  been  obtained 
if  a  larger  number  of  histories  had  been  used  than  the  4000  histories 


which  were  used  in  each  of  the  LITE-II  calculations. 


TRANSMITTED  SCATTERED  INTENSITY  (photons cm_Jsterodion-' /source  photon  cm 
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Fig. IS.  Diffuse  Transmission  and  Reflection  for  Case  E 


REFLECTED  INTENSITY  (photons  cm_Jsterodian_'/source  photon  ci 


Fig,  16.  Diffuse  Tronsmission  ond  Reflection  for  Case  H 
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IV.  CLOUD  ALBEDO  STUDIES 

Monte  Carlo  calculations  were  run  using  the  LITE-II  code  to  deter¬ 
mine  the  number  albedo  for  0.45  micron  light  incident  to  the  top  of  a 
1  km  thick  cloud  at  angles  of  incidence  of  0,  30,  60  and  75  degrees. 

The  phase  function  ':or  0.45  micron  light  scattering  by  the  aerosol 
particles  contained  in  the  cloud  was  taken  from  the  data  reported  in 
Reference  13  for  an  aerosol  size  distribution  that  had  been  suggesced 
by  Deirmendjian  as  typical  of  that  found  in  cumulus  clouds.  The  Monte 
Carlo  results  were  obtained  from  sample  sizes  of  2000  histories  per 
source  angle  with  a  maximum  of  100  collisions  per  history.  The  Monte 
Carlo  calculated  scattered  intensities  in  a  given  polar  angle  interval 
summed  over  all  azimuthal  angles  of  reflection  were  multiplied  by 
[cos0i_^  +  COS0J 

4u(cos0i_^  -  cos0^) 

where  0^  ^  and  0.^  are  the  lower  and  upper  bounds  of  the  polar  angle 
interval  to  give  the  number  of  photons  reflected  per  steradian  in 
direction  3  per  unit  flux  incident  on  the  cloud. 

Figures  17  through  20  present  the  calculated  number  albedo  as  a 

function  of  the  polar  angle  oi  reflection  for  each  angle  of  incidence 

considered.  The  smooth  curves  shown  in  Figures  17  through  20  are 

cosine  distributions  normalized  so  that  the  integral 

*/2 

2irN(0£i,0  *  0°)  /  cos0sin0d0  , 

is  equal  to  the  Monte  Carlo  calculated  value  of  the  total  number  albe¬ 
do,  where  N(0o>0  ■  0°)  is  the  intercept  at  0  ■  0®  for  the  cosine 
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",  POLAR  ANGLE  (deg) 


Fig.  17.  Folor  Angle  Dislribution  of  the  Number  Albedo  for  0.45 p  Light  Incident  on  o  I  Km  Thick 
Cumulus  Cloud:  "Q  0° 
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0,  POLAR  ANGLE  (deg) 


Fig.  18.  Polar  Angle  Distribution  of  the  Number  Albedo  for0.45jx  Light  Incident  on  o  1  Km  Thick 
Cumulus  Cloud:  0O  =  30° 


NUMBER  ALBEDO  (photons  cm"’  ster"  7  photon  cm-3 
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Fig.  19.  Polor  Angulo;  Distribution  of  the  Number  Albedo  for  0.45^  Light  Incident  on  o  1  Km  Thick 
Cumulus  Cl«»d:  tt0  -  60° 


NUMBER  ALBEDO  (photons  cm"1  ster-l/photon 
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Fig.  20.  Polar  Angulor  Distribution  of  the  Number  Albedo  for  0.45/c  Light  Incident  on  a  1  Km  Thick 
Cumulus  Cloud:  l>  75° 

O 
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distribution.  It  is  seen  that  the  angular  distributions  as  obtained  from 
the  Monte  Carlo  calculations  are  closely  approximated  by  cosine  distribu¬ 
tions. 

The  total  number  albedo  as  a  function  of  the  angle  of  incidence  was 
found  to  be  closely  approximated  by  the  expression 

a(0  )  *  0.639  -  O.154cos0 
o  o 

where  0  is  the  angle  of  incidence.  The  results  of  Monte  Carlo  calcula- 
o 

cions  reported  by  R.  E.  Maerker  (Reference  12)  for  the  thermal  neutron 
number  albedo  of  concrete  gave  the  same  functional  dependence  on  the 
angle  of  incidence. 

The  differential  number  albedo  for  angles  of  incidence  other  than 
0q  ■  0°  was  found  to  be  dependent  on  the  azimuthal  angle  of  reflection. 
Figure  21  illustrates  this  dependence  for  0q  =  75°.  The  data  shown  in 
Figure  21  give  the  dependence  of  differential  number  albedo  on  the  polar 
angle  of  reflection  in  the  plane  of  incident 

The  dependence  of  the  intensity  reflected  from  a  cloud  on  the  order 
of  collision  is  shown  in  Figure  22.  The  rate  at  which  the  sum  of  the 
reflected  intensity  increases  with  collision  number  is  seen  to  be  de¬ 
pendent  on  the  angle  of  incidence  for  collision  numbers  less  than  30. 

For  collision  numbers  greater  than  30,  the  curves  have  the  same  shape, 
indicating  that  the  contribution  to  che  sum  of  the  reflected  intensities 
by  a  given  order  of  collision  is  independent  of  the  angle  of  incidence. 
One  can  infer  from  these  data  that  the  fraction  of  the  total  reflected 
intensity  that  has  arisen  from  orders  of  scattering  greater  than  30  is 
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COS  9 


fig.  21 .  Polar  Angular  Distribution  of  the  Differential  Number  Albudo  in  the  Plane  of  Incidence:  90  -  75° 


:e  photon  c 


COLLISION  NUMBER 


Flq.  22.  Reflected  Intensities  from  n  I  Km  Thick  Cumulus  Cloud  Summed  over  Collision  Number 
foi  Vorious  Angles  of  Incidence:  0.45ri  Light 


independent  of  the  shape  of  the  phase  function.  The  leveling  off  of  the 
cumulative  sums  o:  e  reflected  intensities  with  collision  number  indi¬ 

cates  that  contributions  to  the  reflected  intensity  by  orders  of  scatter¬ 
ing  greater  than  100  would  be  negligible.  It  appears  that  less  than  a 
10  percent  error  would  have  been  introduced  in  the  results  if  the  his¬ 
tories  had  been  terminated  after  70  collision,  thereby  reducing  the 
machine  time  requirements  for  the  calculations  by  one  third. 
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V.  SENSITIVITY  OF  LIGHT  SCATTERING  CALCULATIONS 
TO  .CHANGES  IN  ATMOSPHERIC  PARAMETERS 

The  LITE-I  code  was  developed  for  use  in  computing  the  transport  in 
the  atmosphere  of  the  thermal  radiation  emitted  by  a  nuclear  weapon  deton¬ 
ation.  The  wide  range  of  atmospheric  conditions  for  which  transmission 
data  is  needed  indicated  that  a  study  should  be  performed  to  investigate 
the  sensitivity  of  the  LITE-I  results  to  changes  in  the  several  pa.ameters 
that  define  ar.  atmosphere..  The  parameters  selected  for  study  were  the 
source  and  receiver  altitude,  the  aerosol  number  density,  and  the  aerosol 
size  distribution.  Changes  in  the  number  density  affect  only  the  nagni- 
tude  of  the  aerosol  attenuation  coefficient  whereas  changes  in  the  size 
distribution  affect  both  the  magnitude  of  the  aerosol  attenuation  coef¬ 
ficient  and  the  shape  of  the  aerosol  phase  function.  Therefore,  it  is 
useful  to  know  how  much  change  would  occur  in  calculations  of  the  scat¬ 
tered  intensities  at  a  given  receiver  point  if  the  aerosol  number  density, 
size  distribution  and  altitude  variation  were  varied. 

5< 1  Dependence  on  Source  and  Receiver  Altitudes 

Since  the  extinction  cross  section  for  a  given  wavelength  of  light 
varies  with  altitude,  it  is  useful  to  know  how  the  atmospheric  attenuation 
of  light  emitted  by  a  point  isotropic  source  varies  with  both  the  source 
altitude  and  receiver  altitude.  To  investigate  the  dependence  of  the  light 
intensity  on  both  receiver  and  source  altitude,  several  LITE-I  problems 
were  run  for  source  altitudes  of  1,  10,  20,  30  and  45  km.  Point  receivers 
were  placed  at  horizontal  ranges  to  75  km  in  a  horizontal  plane  through 
the  source  and  in  a  horizontal  plane  located  at  an  altitude  of  0. 1  km. 
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The  Elterman  model  atmosphere  was  used  to  describe  the  variation  with 

altitude  and  wavelength  of  the  Rayleigh  scattering,  aerosol  scattering, 

and  ozone  absorption  coefficients.  The  shape  of  the  aerosol  phase  function 

was  obtained  from  an  integration  of  Mie  data  for  spherical  particles  with 

an  index  of  refraction  of  1.5  over  a  particle  size  distribution  given  by 

-4 

the  equation  N(r)  %  r  (Reference  13). 

The  atmospheric  transmission  calculations  were  made  for  monochromatic 
sources  with  wavelengths  of  0.45  and  0.65  microns.  The  LITE-I  output  gave 
the  scattered  Intensities  at  the  receivers  for  ground  albedos  of  0.0  and 
0.9.  The  angular  distribution  of  the  ground  albedo  was  taken  to  be  that 
for  a  Lambert-type  surface. 

The  direct  intensities  at  the  various  receiver  positions  were  computed 
by  use  of  the  equation 

-p (HS,HD,HR) 

I  (HS.HD.HR)  -  - - - - =- 

4ir[(HS-HD)  +HR  ] 

where  HS  is  the  source  altitude, 

HD  is  the  receiver  altitude, 

HR  is  the  horizontal  range,  and 

p(HS,HD,HR)  is  the  mean-free-path  distance  between  the  source  and  receiver. 
The  mean-free-path  distance  p(HS,HD,HR)  was  computed  by  use  of  the  equation 

p (HS,HD,HR)  -  ^  /  p*  (h)dh 

where  p'(h)  is  the  mean-free-path  distance  from  ground  level  to  height  h. 

The  results  obtained  for  the  direct  intensity  at  each  of  the  receiver  points 
are  presented  in  Table  II. 


Table  11=  Direct  Intensities  for  Point  Isotropic  Sources 
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The  LITE-I  calculated  scattered  Intensities  were  added  to  the  direct 
Intensities  to  obtain  the  total  Intensity  at  each  receiver  position.  The 
results  are  presented  In  Table  III  as  a  function  of  the  source  wavelength, 
ground  albedo,  receiver  altitude,  source  altitude  and  horizontal  range. 

2 

The  total  Intensities  given  in  Table  III  were  multiplied  by  4itR  where 
R,  the  source-receiver  distance  in  km,  is  given  by  the  equation 
R  -  [(HS-HD)2  +  HR2]5*  . 

The  results  were  then  plotted  as  a  function  of  R  for  each  source  altitude. 

2 

Multiplying  the  tocal  Intensities  by  4irR  removes  the  geometrical  attenua¬ 
tion  and  reduces  the  numerical  range  over  which  the  intensities  vary. 

2 

The  results  obtained  for  the  0.45  micron  light  sources  giving  4irR  times 
the  total  intensity  as  a  function  of  the  slant  range  R  are  presented  in  Fig¬ 
ures  23  and  24  for  a  receiver  altitude  of  0.1  km  and  in  Figures  25  and  26 
for  the  case  where  the  source  and  receivers  were  at  the  same  altitude.  Simi¬ 
lar  data  for  the  0.65  micron  light  sources  are  presented  in  Figures  27  through 
30. 

A  comparison  of  the  data  shown  in  Figures  23  through  30  for  ground  albe¬ 
dos  of  0.0  and  0.9  show  that  the  effect  of  increasing  the  ground  albedo  is 
to  increase  the  total  intensity  at  a  given  receiver  position  over  that  com¬ 
puted  for  a  ground  albedo  of  0.0.  The  amount  of  the  increase  in  the  inten¬ 
sity  resulting  from  ground  scattering  is,  for  a  given  value  of  the  ground 
albedo,  a  function  of  the  source  altitude,  the  receiver  altitude,  and  the 
horizontal  range.  The  dependence  of  the  intensity  at  a  given  receiver  posi¬ 
tion  on  the  magnitude  of  the  ground  albedo  was  found  to  be  exponential. 
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R,  SUNT  RANGE  (km) 


Fig.  23.  Variotion  of  An  RJ  Times  Totcl  Intensity  with  Slont  Konge:  0.45/t  Light,  Ground  Albedo  0.0, 
Receiver  Altitude  ~  0. 1  Km 


4ir  RJ  TIMES  TOTAL  INTENSITY  (photons/source  photon) 
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R,  SLANT  RANGE  (km) 


Fig.  24.  Vcriatlon  of  4wRs  Times  Total  Intensity  with  Slant  Range:  0,45] u  Light,  Ground  Albedo  =  0.9, 
Receiver  Altitude  =  0.1  Km 


4trR3  TIMES  TOTAL  INTENSITY  (photons/source  photon) 
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Fig.  25.  Variation  af  4irR3  Times  Total  Intensit  with  Slant  Range:  0.45/*  Light,  Receiver  and  Source 
at  Same  Altitude,  Ground  Albedo  =  0.0 


47tR2  TIMES  TOTAL  INTENSITY  (phofons/soorce  photon) 
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R,  SLANT  RANGE  (km) 


Fig.  26.  Vorioticn  of  4wRJ  Times  Totol  Intensity  with  Slont  Range:  0 . 45m-  Light,  Receiver  and  Source 
at  Same  Altitude,  Ground  Albedo  =  0.9 


4rr  R3  TIMES  TOTAL  INTENSITY  (photon s/sourc.3  photon) 
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R,  SLANT  RANGE  (km) 


rig.  27.  Variotion  of  4w  R3  Times  Total  Intensity  with  Slont  Ronge:  0.65/1  Light,  Receiver  Altitude  =  0.1  Km, 
Ground  Albedo  =  0.0 


Air  R2  TIMES  TOTAL  INTENSITY  (photons/source  photon) 
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R,  SLANT  RANGE  (km) 


Fig.  28.  Variation  of  Air  R*  Times  Total  Intensity  with  Slant  Range:  0.65 (i  Light,  Receiver  Altitude  =0.1  Km, 
Ground  Albedo  =  0.9 


4irR/  TIMES  TOTAL  INTENSITY  (photons./soorce  photon) 
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Fig.  29.  Variation  of  4wR7  Times  Total  Intensity  with  Slant  Range:  0.65 /t  Light,  Receiver  and  Source 
at  Same  Altitude,  Ground  Albedo  =  0.0 


4-77- R 2  TIMES  TOTAL  INTENSITY  (photons/source  photon) 
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0  10  20  30  40  50  60  70  80  . 

R,  SLANT  RANGE  (K<") 


Fig.  30.  Variation  of  4tR2  Times  Total  Intensity  with  Slant  Range:  0.65/4  Light,  Receiver  and  Source 
at  Same  Altitude,  Ground  Albedo  =  0.9 
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2 

It  is  seen  in  Figures  23,  24,  27,  and  28  that  4nR  times  the  total  in¬ 
tensity  increases  with  the  increasing  source  altitude  when  the  receivers  are 

located  near  the  ground  surface.  This  increase  in  the  total  intensity  times 
2 

4nR  results  mainly  from  the  fact  that  the  mean-free-path  distance  for  a 
given  slant  range  decreases  as  the  source  altitude  is  increased. 

2 

The  variation  of  4nR  times  the  total  intensity  with  altitude  for  ground 
ranges  of  0,  10,  20,  50,  and  75  km  is  shown  in  Figure  31  for  0.45  micron 
light  and  a  ground  albedo  of  0.0.  For  a  ground  range  of  0.0  the  total  inten¬ 
sity  times  4itR^  decreases  as  the  source  altitude  increases.  For  ground 

2 

ranges  of  20  km  and  greater,  the  total  intensity  times  4ttR'  increases  as 

the  source  altitude  increases.  The  rate  at  which  the  total  intensity  times 
2 

4irR  increases  with  source  altitude  for  horizontal  ranges  greater  than  10  km 

is  seen  in  Figure  31  to  be  dependent  on  the  horizontal  range.  The  magnitude 

2 

of  the  total  intensity  times  4ttR  for  receivers  placed  near  the  ground  and 
at  ground  ranges  up  to  75  km  appears  to  be  approaching  a  value  which  is 
independent  of  ground  range  as  the  source  altitude  is  increased. 

The  variation  of  the  total  intensity  with  horizontal  range  for  a  source 
wave  length  of  0  45  microns  and  a  ground  albedo  of  0.0  is  shown  in  Figure  32 
for  each  source  altitude  studied.  For  a  given  horizontal  range  the  maximum 
intensity  is  seen  to  be  dependent  on  the  source  altitude.  The  altitude  at 
which  the  maximum  intensity  is  obtained  is  seen  in  Figure  32  to  increase 
with  increasing  horizontal  range  for  ranges  less  than  80  km. 

The  meah-f ree-path  distance  for  a  given  source-receiver  distance  de¬ 
creases  as  the  receiver  altitude  is  raised  to  the  source  altitude.  This  effect 


4ttRj  TIMES  TOTAL  INTENSITY  (photons/source  photon) 
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Fig.  31,  Variation  of  4  7tR 2  Times  Total  Intensity  with  Source  Altitude:  0.45/u  Light,  Receiver  Altitude  - 
0.1  Km,  Ground  Albedo  =0.0 


INTENSITY  (photons  km"J/source  photon) 
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Fig.  32.  Total  Intensity  vs  Horizontal  Range  for  Several  Source  Altitudes:  0.45m  Light, 
Receiver  Altitude  =  O.lKm,  Ground  Albedo  =  0.0 
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is  seen  in  Figure  33  where  4itR  times  the  total  intensity  as  a  function 
of  the  source- receiver  distance  is  shown  for  receivers  located  at  an  alti¬ 
tude  of  0.1  km  and  at  the  same  altitude  as  the  source.  The  data  shown  in 
Figure  33  are  for  0.45  y  light,  a  ground  albedo  of  0.0,  and  source  alti¬ 
tudes  of  10  and  30  km.  It  is  seen  that  for  a  given  slant  range  and  source 
altitude,  the  intensity  increases  as  the  receiver  altitude  is  increased  up 
to  the  source  altitude. 

5.2  Dependence  on  Aerosol  Number  Density 

The  effect  of  varying  the  aerosol  number  density  on  the  scattering  of 

light  in  the  atmosphere  was  studied  for  atmospheres  defined  by  applying 

perturbations  to  an  atmosphere  for  which  only  Rayleigh  scattering  and  ozone 

absorption  were  allowed.  The  cross  sections  for  Rayleigh  scattering  and 

ozone  absorption  were  taken  from  a  tabulation  reported  reported  by  Elterman 

(Reference  2) .  This  basic  atmosphere  was  then  perturbed  by  the  addition 

of  aerosol  particles  with  number  densities  of  2300,  4600,  and  9200  part- 
3 

icles  per  cm  at  sea  level  air  density  and  with  a  size  distribution  defined 
by  Deirmendjian' s  (Reference  14)  Haze  C  model  for  continental  aerosols. 

The  vertical  distribution  of  the  aerosol  particle  concentration  was  assumed 
to  be  the  same  as  that  formulated  by  Elterman  from  field  measurements  re¬ 
ported  by  Penndorf  (Reference  15)  and  by  Chagnon  and  Jun^e  (Reference  16). 

The  aerosol  cross  section  at  sea  level  air  density  for  Haze  C  type  aero¬ 
sols  was  obtained  by  integrating  Mie  scattering  data  for  spherical  part¬ 
icles  with  an  index  of  refraction  of  1.5  over  the  Haze  C  size  distribution 
(Reference  13).  The  variation  of  the  ground  level  aerosol  scattering  co¬ 
efficient  for  the  Haze  C  size  distribution  with  wavelength  is  given  in 
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Fig.  33.  Dependence  of  4ir  Times  Totol  Intensity  on  Receiver  Altitude:  0.45/r  Light,  Ground 
Albedo  0.0 
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Table  IV  for  a  number  density  of  one  particle  cm  .  Aerosol  scattering 
coefficients  for  a  given  number  density  of  Haze  C  aerosols  can  be  obtained 
by  multiplying  the  data  in  Table  IV  by  the  aerosol  number  density  in  part¬ 
icles  cm  \ 


Table  IV.  Scattering  Coefficient  for  Haze  C  Aerosols 
With  a  Number  Density  of  One  Particle  cm--* 

-1  -3 

(km  /particles  cm  ) 


Wavelength 

(y) 

Aerosol  Scattering 
Coefficient 

0.3 

1.090-04 

0.45 

7.239-05 

0.5 

6.574-05 

0.65 

5.200-05 

0.7 

4.891-05 

Normalized  aerosol  scattering  phase  functions  for  several  light  wave¬ 
lengths  were  also  computed  for  the  Haze  C  size  distribution.  The  results 
are  given  in  Table  V  as  a  function  of  the  scattering  angle. 

All  of  the  calculations  performed  for  use  in  this  study  were  for  point 
isotropic  monochromatic  sources  located  at  an  altitude  of  807.7  m.  The 
receivers  were  placed  at  an  altitude  of  6.1  m  and  at  horizontal  ranges  of 
up  to  48.3  km.  The  source  wavelengths  studied  were  0.3,  0.45,  0.5,  0.65, 
and  0.7p. 

The  printed  output  of  the  LITE-I  calculations  gave  the  scattered  and 
direct  intensities  at  a  spherical  detector.  The  detector  always  presented 
a  unit  area  in  a  direction  normal  to  the  direction  of  motion  of  the  photons 


xzmamm 
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TABLEV.  NORMAL  1 7  FO  VOLUMfc  SCATTERING  FUNCTION  :  HA7F  C  NOOFL 
INDEX  OF  REFRACT ION  *  l.SO 


SCATTERING 


ANGLE 

WAVELENGTH  (M1CKUNS) 

(deg) 

.  10 

.45 

.50 

.65 

.70 

Oc 

7.  2  7  OF  CO 

4.627F  00 

4.161F  00 

3 . 8 1 1 F  00 

3.493F  00 

1.0 

5.606F  00 

4. 089F  00 

3.736F  00 

1.513F  00 

3.261F  00 

7.0 

3.468F  00 

3.034E  00 

7.859F  00 

7.837F  00 

2.713F  00 

3.0 

2.575E  00 

7.74  7F  00 

2. 14  3F  00 

7.183F  00 

2.142E  00 

4.0 

7.016F  00 

1.815F  00 

1.777F  00 

1.749F  00 

1.728F  00 

5.0 

1.696F  00 

1.5  3  3F  00 

1.464F  00 

1.485E  00 

1.466F  00 

6.0 

1.469F  00 

1.376F  00 

1.770F  00 

1.798F  00 

1.283F  00 

7.0 

1.798F  00 

1.1 72F  00 

1.I75F  00 

1.151F  00 

1.138F  00 

8.0 

1.163F  00 

1.05  IF  00 

1.009F  00 

1.036F  00 

1.077E  00 

9.0 

1.059F  no 

9.518F-01 

9.  164F-01 

9.450F-01 

9. 332F-01 

10.0 

9.  744F-01 

8. 757F-01 

8.477F-01 

8.731E-01 

8.638F-01 

1  1.0 

9.050F-01 

8. 149F-01 

7.86 1F-01 

8.  167E-01 

8.087E-01 

12.0 

8.519F-01 

7. 6  78F-01 

7.41 7F-01 

7.688F-01 

7. 604F-01 

13.0 

7.994F-01 

7. 722F-01 

6. 978F-01 

7.212F-01 

7. 135E-01 

14.0 

7.  380F-01 

6. 733F-01 

6. 5 1 6F-0 1 

6.7  1 5F-01 

6.648F-01 

15.0 

6.870F-01 

6. 761F-01 

6.070F-01 

6. 2  33F-01 

6. 175F-01 

16.0 

6. 367F-01 

5. 866F-0 1 

5.690F-01 

5. 819F-01 

5. 765F-0 1 

1  7.0 

5.970F-01 

5. 537F-01 

5. 373F-01 

5.486E-01 

5.432F-01 

18.0 

5.640F-01 

5. 777F-01 

5.087F— 01 

5. 204F-01 

5. 157F-0 1 

19.0 

5.336F-01 

4.954E-01 

4.875F-01 

4.940E-01 

4. 895E-01 

20.0* 

5.  051F-01 

4.704F-01 

4.586F-01 

4.687E-01 

4.647F-01 

27.5 

4.  387F-01 

4.10  7F-01 

4.01 9F-01 

4.  1  11F-01 

4.075F-01 

25.0 

3. 781 F -01 

3.  586F-01 

3.52  3F-01 

3. 581F-01 

3. 554F-0 1 

27.5 

3.769F-01 

3. 156F-01 

3.  1 1  OF -0 1 

3.  129E-01 

3. 108F-01 

30.0 

7.85  7F-01 

7.796F-01 

7. 764F-01 

7. 763F-01 

7.744F-01 

37.5 

2.479F-01 

2. 468E-01 

7.450F-01 

7.428F-0 1 

2.416E-01 

35.0 

2.177F-01 

7. 1 8 1F-01 

7. 174F-01 

?.  1 12F-01 

2. 102F-01 

3  7.5 

1.8S3F-01 

1 . 929F -Oi 

1.931F— 01 

1.865F-01 

1.858E-01 

40.0 

1.678F-01 

1.706F-01 

1.715F-01 

1.658E-01 

1 . 655F-0 1 

42.5 

1 . 409F-01 

1.51 7E-01 

1.525F-01 

1 . 458F-01 

1.455F-01 

45.0 

1.2  32F-01 

1. 336F-01 

1. 354F-01 

' • 2  93F-0 1 

1.292F-01 

4  7.5 

1.090F-01 

1.  182F-01 

1. 704F-01 

1.158F-01 

1. 159F-01 

50.0 

9.577F-0? 

1.054F-01 

1 . 078F-01 

1.078E-01 

1.030E-01 

57.5 

8.3  36F-0? 

9.  324F-07 

9.  580F-07 

9.  1  18F-02 

9.  162F-07 

55.0 

7. 36  IF- 07 

8.  766F-02 

8.  5 3 IF  -07 

8. 1 78F-02 

8.714E-02 

5  7.5 

6.  567F-02 

7.451F-02 

7.716F-02 

7.375F-02 

7.41 3E-0? 

60.0 

5.  704F-0? 

6.  589F-02 

6. 855F-02 

6. 570F-02 

6. 564F-07 

67.5 

5.  1  1  IF -07 

5. 974F-02 

6.  183F-07 

5. 900F-07 

5.951F-02 

65.0 

4.561F-0? 

5.  306F-02 

5.558F-02 

5. 3  30F-02 

5. 382E-0? 

67.5 

4.  063F  - 07 

4. 7  77F-07 

5.01 6F-07 

4.806F-02 

4. 858F-02 

70.0 

3.675F-0? 

4.307F-07 

4.54IF-02 

4, 3  37F-02 

4. 397E-07 

77.5 

3.  764F-07 

3.873F-07 

4.097F-07 

3 . 941 E -02 

3. 998E-07 

75.0 

7.966F-0? 

3.576F-02 

1.718F-07 

3.605F-02 

3.659F-02 

7  7.5 

7.658F-07 

1. 197F-0? 

3. 397F-07 

3.767F-0 2 

3. 327F-07 

80.0 

2.4  30F -07 

2.91 1F-07 

3.099F-07 

3 • 003F-02 

3.059F-07 

87.5 

7.740F-02 

7. 6  77F-02 

7.848F-02 

2. 774F-02 

?• 879F-07 

85.0 

7.046F-07 

7.444F-07 

7.61 0F-02 

2.548F-07 

2  603F-02 

87.5 

1.858F-07 

7.221F-02 

2. 179F-07 

7. 334F-02 

7. 390F-02 

90.0 

1. 687F-07 

7.01  7F-07 

7. 166F-07 

2.142F-02 

2. 199F-07 

97.5 

1. 555F-0? 

1. 864F-02 

7.003F-0? 

1 . 987F-02 

7.045F-07 

69 


TABLE  V.  NORMAL  12 ED  VOLUME  SCATTERING  EUNCTION  :  HAZE  C  MODEL 
INUEX  Of  kEFRACT l UN  >  1.50 


SCATTERING 

ANGLE 

(deg)  .  30 

95.0  1.466E-02 

97.5  J.409F-02 

100.0  1.328E-02 

102.5  1.255F-02 

105.  0  1.  1866-02 

107.5  1. 112E-02 

110.0  1.060E-02 

112.5  1.038F— 02 

115.0  1.024F-02 

117.5  1.010E-02 

120.0  I.017E-02 

122.5  1.025E-02 

125.0  1.027F-02 

127.5  1.045F-02 

130.0  1.062E-02 

132.5  1.086E-02 

135.0  1.U5E-02 

137.5  1.150E-02 

140.0  1.215E-02 

142.5  1.292E-02 

145.0  I.373E-02 

147.5  1.502E-02 

150.0  1.653E-02 

152.5  1. 807E-02 

155.0  2.035F-02 

156.0  2.123E-02 

157.0  2.206E-02 

158.0  2.303E-02 

159.0  2.392E-02 

160.0  2.460E-02 

161.0  2.493F-02 

162.0  2.540E-02 

163.0  2.616E-02 

164.0  2.648E-02 

165.0  2.6r  8E-02 

166.0  2.663E-0? 

167.0  2. 714F-02 

168.0  2. 785E-02 

169.0  2.840F-02 

170.0  2.897E-02 

171.0  2.965E-02 

1/2.0  3.060E-02 

173.0  3. 164F-02 

174.0  3. 304F-02 

175.0  3.460E-02 

176.0  3. 610F-02 

177.0  3.683F-02 

178.0  3. 721E-02 

179.0  3. 521E-02 

180.0  3.534F-02 


WAVELENGTH  {MICRONS) 

.45  .50 

1 • 74  If -02  1.871F-02 

l. 6536-02  1. 774E-02 

I.55IF-02  1.665E-02 

1.455F-02  1.563F-02 

t. 368E-02  1.471F-02 

1. 274E-02  1.374E-02 

1.21 3E-02  1.308E-02 

1.172E-02  1.26  3E-02 

I. 143E-02  1.229E-02 

1.  1246-02  1.207E-02 

1.  1 14E-02  I. 1946-02 

1. 108F-02  1.184E-02 

1. 1036-02  1. 176F-02 

1.  103E-02  l. 175E-02 

1. 109E— 02  1.178E-02 

1. 1286-02  1 . 195E-02 

1. 163E-02  1.227E-02 

1. 207F-02  1. 268E-02 

1.265E-02  1.322E-02 

1. 3296-02  1. 3836-02 

1.395E-02  1.445E-02 

1. 494E-07  1.539E-07 

!. 6176-02  1.655E-02 

1. 755E-02  I. 785F-02 

1.941E-02  1.962F-02 

2. 020E-02  2.035E-02 

2. 093E-02  2. 105E-02 

2.169E-02  2.173E-02 

2.242E-02  2.24  IE-02 

2.3I3F-02  2.312E-02 

2.381E-02  2.374E-02 

2.427E-02  2.418E-02 

2.4566-02  2.447E-0 2 

2.486E-02  2.473E-02 

2.506E-02  7.496E-02 

2.526E-02  2.516E-02 

2.565E-02  2.551E-0? 

2.612E-02  2.598E-02 

2. 656E-02  2.636E-02 

2. 701F-02  2.680E-02 

2. 776E-02  2.755E-02 

2. 87BF-02  2.844E-02 

2.98IE-02  2.948F-02 

3. 106F-02  3.062F-02 

3.219F-02  3. 163E-02 

3. 334F-02  3. 284E-02 

3.408E-02  3. 332F-02 

3. 318E-02  3. 24 IE -02 

3. 265E-02  3.213F-02 

3. 304E-02  3. 253E-02 


.65  .70 

I.873E-02  1.931F-0? 
1.791E-02  1.849E-02 
1.692E-02  1  *  750E-02 
1.603E-02  1.66  IE-02 
1.521E-02  1.581E-02 
1.437E-02  1.499E-02 
1. 376E-02  1.438E-02 
1 . 344E-02  1.407E-02 
1. 321E-02  1. 385E-02 
1.299E-02  1.36E-02 
1.298E-02  1.364E-02 
1. 300E-02  1.367E-02 
1 . 298E-02  1.367E-02 
1. 310E-02  1. 380E-02 
1. 325E-02  1.396E-02 
1.345E-02  1.417E-02 
1.370E-02  1. 444E-02 
U402E-07  1.476E-02 
l .460E-02  1. 535E-02 
1.528E-02  1.605F-02 
* . 601 E-02  1.678F-02 
1. 715E-02  1. 792E-02 
I.849E-02  1.924E-02 
l. 982E-02  2.055E-02 
2.171F-02  2.243E-02 
2.242E-02  2.312E-02 
2.305E-02  2.375E-02 
2.378E-02  2.445E-02 
2.431E-02  2.513E-02 
2.497E-07  2.5606-02 
2.525E-02  Z.593E-02 
2. 574F-02  2.639E-02 
2.643E-02  2.703E-02 
2.690E-02  2.753E-02 
2.702F-02  2.767E-02 
2.712E-02  2.776F-02 
2.751E-02  2.816E-02 
2.804E-02  2.8746-02 
2.850E-02  2.9176-02 
2.8986-02  2.9546-02 
2.956E-02  3.0166-02 
3.0  32E-02  3.0976-02 
3.126F-02  3.1786-02 
3.2286-02  3.2846-02 
3. 362F-02  3. 4296-02 
3.486F-02  3. 526F-02 
3.4816-02  3.4976-02 
3. 398E-02  3.4386-02 
3.4 10F-02  3.4706-02 
3.455F-02  3.513E-02 
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entering  the  detector,  irrespective  of  their  direction  of  motion.  The 
scattered  intensities  were  also  printed  out  as  a  function  of  the  ground 
albedo  for  ground  albedos  of  0.0  and  C.9.  The  ground  was  treated  as  a 
Lambert  type  reflection  surface. 

Table  VI  presents  the  calculated  direct  intensities  at  each  of  five 
different  receiver  positions  as  a  function  of  the  source  wave  length  and 
the  aerosol  number  density.  Table  VII  presents  the  calculated  scattered 
intensities  at  the  five  different  receiver  positions  for  ground  albedos 
of  0.0  and  0.9  as  a  function  of  the  source  wave  length  and  the  aerosol 
number  density.  From  an  examination  of  the  scattered  intensity  data  in 
Table  VII,  it  is  seen  that  the  scattered  intensity  at  4.03  km  ground  range 
increases  with  increasing  aerosol  number  density,  whereas  at  larger  ground 
ranges  he  scattered  intensity  decreases  with  increasing  aerosol  number 
density. 

The  variation  of  the  scattered  intensities  at  each  of  five  different 
receiver  positions  with  source  wave  length  is  shown  in  Figure  34  for  the 
basic  atmosphere  (Rayleigh  scattering  plus  ozone  absorption).  The  data 
shown  in  Figure  34  are  for  a  ground  albedo  of  0.9.  Also  shown  in  Figure 
34  is  a  plot  of  the  Rayleigh  scattering  coefficient  at  sea  level  altitude 
versus  wave  length.  In  addition  to  the  source  wave  lengths  listed  in 
Table  VI,  source  wave  lengths  of  U.34  and  0.38u  were  also  used  in  calcula¬ 
tions  to  study  the  effect  of  the  ozone  absorption  on  the  scattered  inten¬ 
sities.  It  is  seen  that  the  strong  ozone  absorption  occurring  near  0.3p 
results  in  a  rapid  decrease  in  the  scattered  intensities  as  the  source 
wave  length  approaches  0.3u  for  ground  ranges  greater  than  4.03  km.  For 
wave  lengths  greater  than  0.34y  the  scattered  intensities  at  each  receiver 
varies  approximately  the  same  with  wave  length  as  does  the  Rayleigh 
scattering  coefficient. 


Table  VI.  Direct  Intensities  as  a  Function  of  the  Number 
Density  of  Haze  C  Aerosols  at  Ground  Level 
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Fig.  34.  Scattered  Intensities  vs  Source  Wavelength  far  Several  Ground  Ranges: 
Rayleigh  Atmosphere  with  Ozone  Absorption 
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2 

Figure  35  presents  the  variation  in  the  scattered  intensity  time  4itR 
at  several  ground  ranges  with  the  aerosol  cross  section  at  sea  level  alti¬ 
tude  for  a  source  wavelength  of  0.45y  and  a  ground  albedo  of  0.9.  It  is 
seen  that  at  close-in  distances  the  scattered  intensity  for  a  point  iso¬ 
tropic  0.45y  light  source  increases  with  increasing  aerosol  cross  section. 

At  larger  ground  ranges  the  scattered  Intensities  decrease  when  the  aerosol 

2 

cross  section  is  increased.  The  scattered  intensity  time  4uR  for  an  aero¬ 
sol  cross  section  of  0.180  km  ^  corresponds  to  the  attenuation  that  would 
be  expected  for  an  atmosphere  defined  by  the  "clear  standard  atmosphere" 
formulated  by  Elterman. 

Figure  35  also  shows  the  variation  with  the  aerosol  cross  section  of 

2 

the  total  intensity  (scattered  plus  direct  intensity)  time  4itR  for  0.45y 
light  at  ground  ranges  to  48.38  km.  The  ground  albedo  was  0.9.  Except  at 
ground  ranges  less  than  12.07  km,  the  total  intensity  results  almost  en¬ 
tirely  from  scattered  radiation  when  the  aerosol  cross  section  is  greater 

than  0.2  km  It  is  seen  from  the  data  in  Figure  35  that  the  variation 

2 

of  the  total  intensity  time  4ttR  with  ground  range  is  highly  dependent  on 
the  magnitude  of  the  aerosol  cross  section  for  values  of  the  aerosol  cross 
section  up  to  twice  the  value  used  by  F.lterraan  for  a  "clear  standard  atmos¬ 
phere."  For  larger  values  of  the  aerosol  cross  section,  the  total  inten¬ 
sity  times  4irR“  at  a  given  ground  range  frox  a  0.45y  point  source  does 
not  vary  significantly  as  the  aerosol  cross  section  is  increased. 

2 

The  dependence  of  both  the  scattered  and  total  Intensities  times  4itR 
on  the  magnitude  of  the  aerosol  cross  section  is  shown  in  Figure  36  for  a 
source  wavelength  of  0.7  and  a  ground  albedo  of  0,0.  It  is  seen  that  the 
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g.  35.  Variation  of  47tR5  Times  Scattered  and  Total  Intensities  with 

Aerosol  Cross  Section  for  Several  Ground  Ranges:  0.45/t  Light 


4  7rR2  TIMES  INTENSITY  (m2  photons  m“2/source  photon) 
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CROSS  SECTION  (km'1) 


Fig.  36,  Variation  of  4itR2  Times  Scattered  and  Total  Intensities  with  Aerosol 
Cross  Section  for  Several  Ground  Ranges:  0.70fi  Light 
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dependence  on  the  aerosol  content  of  the  atmosphere  fcr  0.7p  light  is  simi¬ 
lar  to  that  shown  in  Figure  35  for  0.45u  light.  Similarly  it  was  found  for 
the  other  source  wavelengths  studied  (0.3,  0.5,  and  0.65  )  that  the  attenua¬ 
tion  of  both  the  total  and  scattered  Intensities  is  highly  dependent  on  the 
aerosol  number  density.  The  total  intensity  was  found  to  vary  by  as  much 
as  a  factor  of  30  for  a  ground  range  of  48  km  as  the  aerosol  number  density 
is  varied  from  0  to  9200  particles  per  cc. 

The  total  intensity  at  a  given  receiver  position  is  calculated  by  the 
Monte  Carlo  code  as  the  sum  of  the  intensities  received  after  each  colli¬ 
sion.  The  number  of  collisions  that  have  to  be  followed  for  a  given  wave¬ 
length  before  the  contribution  from  high  order  collisions  is  negligible 
is  a  function  of  both  the  aerosol  cross  section  and  the  aerosol  phase  func¬ 
tion  used  in  the  calculation.  Figure  37  shows  for  a  ground  albedo  of  0.9 
the  sum  of  the  scattered  intensities  vs  collision  number  at  a  ground  range 
of  32.19  km  for  0.3p  wavelength  light.  An  increase  in  the  aerosol  cross 
section  at  sea  level  results  in  an  increase  in  the  number  of  collisions 

o 

that  has  to  be  followed  before  the  contribution  to  the  total  scattered 
intensity  by  higher  order  collisions  becomes  negligible.  The  minimum  num¬ 
ber  of  collisions  that  need  to  be  followed  for  an  aerosol  cross  section 
of  0.2506  km  ^  is  approximately  30,  whereas  the  minimum  number  needed  for 
an  aerosol  cross  section  of  0.5012  km  ^  is  40. 

Curves  similar  to  those  shown  in  Figure  37  were  obtained  for  the  other 
source  wavelength  studies.  They  showed  that  an  increase  in  the  sum  of  the 
aerosol  and  Rayleigh  scattering  cross  section  results  in  an  increase  in 


SCATTERED  INTENSITY  (photons  m-2/source  photon) 
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Fig.  37.  Sum  of  the  Scattered  Intensities  vs  Collision  Number  for 
Several  Aerosol  Cross  Sections:  0,30n  Light 


80 


the  number  of  collisions  required  before  the  sum  of  the  scattered  intensi¬ 
ties  resulting  from  higher  order  collisions  becomes  negligible. 

The  ratio  of  the  scattered  intensities  for  the  case  where  the  ground 
albedo  was  0.9  to  that  for  a  ground  albedo  of  0.0  was  found  to  not  vary 
significantly  with  the  horizontal  range  for  all  of  the  source  wavelengths 
studied.  It  was  found  that  the  ratio  decreased  with  increasing  wavelength 
from  the  value  of  about  2.5  for  0.3  microns  to  a  value  of  1.5  for  0.7  microns. 
The  scattered  intensity  for  a  given  horizontal  range  was  found  to  increase 
exponentially  with  an  increase  in  the  ground  albedo. 

5 . 3  Dependence  on  Changes  in  Aerosol  Phase  Function 

To  study  the  effect  of  varying  the  shape  of  the  aerosol  phase  func¬ 
tion  on  the  magnitude  of  the  scattered  light  intensity  in  the  atmosphere, 
LITE-I  problems  were  run  using  various  aerosol  phase  functions.  The  atmos¬ 
phere  used  in  these  problems  was  taken  to  be  that  described  by  Elterman's 
"clear  standard  atmosphere"  model.  The  phase  functions  used  in  the  LITE-I 
calculations  for  aerosol  scattering  were  computed  by  integrating  Mie  scat¬ 
tering  data  for  spherical  particles  with  an  index  of  refraction  of  1,5  with 
particle  size  distributions  defined  by  Diermendjian's  Haze  C  and  Haze  M 
models  (Reference  14)  for  continental  and  maritime  hazes  and  with  exponen¬ 
tial  size  distributions  described  by  N(r)*r  V  for  v  -  2  and  3.  In  some  of 
the  LITE-I  calculations  the  Rayleigh  phase  function  was  used  instead  of  an 
aerosol  phase  function  to  describe  aerosol  scattering. 

The  shape  of  the  normalized  phase  function  for  light  scattering  by 
aerosol  particles  is  dependent  on  the  wavelength  of  the  incident  light. 
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the  particle  size  distribution,  and  the  index  of  refraction.  Phase  func¬ 
tions  for  aerosol  scattering  are  highly  peaked  in  the  forward  direction. 
Table  VIII  lists,  for  wavelengths  of  0.3  and  0.7p  and  for  several  aerosol 
size  distributions,  the  cosine  of  the  upper  bound  of  the  scattering  angle 
interval  into  which  a  given  fraction  of  the  photons  undergoing  collisions 
in  a  unit  volume  of  aerosols  will  be  scattered.  It  is  seen  from  the  data 
in  Table  VIII  that  at  least  90  percent  of  the  photons  undergoing  aerosol 
particle  scattering  in  a  unit  volume  of  air  will  be  scattered  into  angles 
between  0°  and  90°  for  all  of  the  aerosol  size  distribution  models  listed. 

In  contrast,  only  50  percent  of  the  radiation  undergoing  Rayleigh  scat¬ 
tering  will  be  scattered  in  a  forward  direction. 

An  examination  of  the  variation  with  both  wavelength  and  altitude  of 
the  ratio  of  the  Rayleigh  scattering  coefficient  to  the  sum  of  the  Rayleigh 
and  aerosol  scattering  coefficients  for  the  Elterman  "clear  standard  atmos¬ 
phere"  reveals  that  the  probability  of  a  Rayleigh  scattering  event  occur¬ 
ring  at  a  given  altitude  reduces  with  source  wavelength.  It  was  also  found, 
for  wavelengths  of  X  -  0.3,  0.45,  and  0.5p,  that  more  than  90  percent  of 
the  scattering  at  altitudes  above  5  km  are  Rayleigh  events.  For  wavelengths 
of  X  ■  0,65  and  0.70p,  the  probability  of  a  Rayleigh  scattering  event  is 
greater  than  0.9  only  for  altitudes  between  7  and  18  km  and  above  30  km. 
Therefore,  the  use  of  the  Rayleigh  phase  function  to  define  the  phase  func¬ 
tion  for  aerosol  particle  scattering  should  not  have  a  significant  effect 
on  the  calculated  intensities  for  wavelengths  between  0.3  and  0.5  .  For 
wavelengths  between  0.65  and  0.7y,  there  should  be  a  noticeable  difference 
between  a  calculation  which  used  the  Rayleigh  phase  function  to  describe 
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aerosol  scattering  and  a  calculation  which  used  an  aerosol  phase  function 
to  define  the  angular  distribution  of  scattering  by  aerosol  particles. 

LITE-I  calculations  for  use  in  studying  the  effect  of  varying  the 
shape  of  the  aerosol  phase  function  were  run  for  a  point  isotropic  mono¬ 
chromatic  source  located  at  an  altitude  of  807.7  m.  The  receivers  were 
placed  at  an  altitude  of  6.1  m  and  at  horizontal  ranges  of  up  to  48.3  km. 
The  ground  albedo  was  taken  to  be  a  Lambert  type  albedo  with  a  value  of 
0.9.  Results  were  also  obtained  for  a  ground  albedo  of  0.0.  The  source 
wavelengths  studied  were  0.3,  0.45,  0.5,  0.65,  and  0.7p. 

The  calculated  scattered  intensities  at  each  receiver  point  are  pre¬ 
sented  in  Table  IX  as  a  function  of  the  source  wavelength,  the  ground  al¬ 
bedo,  and  the  phase  function  used  to  describe  aerosol  scattering.  The 
direct  intensities  at  each  receiver  point  are  presented  in  Table  X  as  a 
function  of  the  source  wavelength. 

The  results  obtained  from  LITE-I  calculations  for  problems  using 
phase  functions  defined  by  the  Haze  C  and  Haze  M  size  distribution  models 
and  n(r)%r  V  for  v  ■  2  and  3  and  the  Rayleigh  phase  function  for  aerosol 
scattering  are  presented  in  Figures  38  and  39  for  wavelengths  of  0.45  and 
0.7p,  respectively.  For  a  source  wavelength  of  0.45p,  it  is  seen  that  the 
maximum  intensity  computed  for  a  given  horizontal  range  does  not  differ 
by  more  than  35  percent  of  the  minimum  value.  An  examination  of  the  data 
for  0.7p  shows  that  there  are  no  large  differences  between  the  computed 
attenuation  cut  to  50  km  as  given  by  each  of  the  aerosol  models,  but  the 
use  of  the  Rayleigh  phase  function  to  describe  the  aerosol  phase  function 


Table  IX,  Scattered  Intensities  in  an  Eltennan  Atmosphere 

for  Various  Aerosol  Phase  Functions  and  Wavelengths 
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Tabei  X,  Direct  Intensities  in  an  Elterman  Atmosphere  for  Several  Wavelengths 


RJ  TIMES  INTENSITY  (photons /soot ce  photon) 
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ig.  36.  Variotion  of  An  R1  Times  the  Direct  and  Totol  Intensities  with  Horizontal  Range  for  Voriovs 
Aerosol  Scattering  Models:  0.45a?  Light 


4ffRJ  TIMES  INTENSITY  (phctons/source  photon) 
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Fig.  39.  Variation  of  4TR;Times  the  Direct  ondTotol  Intensities  with  Horizontal  Ronge  for  Various 
Aerosol  Scottering  Models:  0 . 7 0/ 1  Light 
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will  give  results  that  are  as  much  as  a  factor  of  two  less  than  that  ob¬ 
tained  using  the  phase  functions  from  the  various  aerosol  models. 

The  scattered  intensities  at  a  given  receiver  position  as  a  function 
of  collision  number  are  dependent  on  the  magnitude  of  the  sum  of  the  aero¬ 
sol  and  Rayleigh  scattering  cross  sections,  the  magnitude  of  the  aerosol 
cross  section,  and  the  shape  of  the  aerosol  phase  function.  The  sum  of 
the  aerosol  and  Rayleigh  scattering  cross  section  for  a  given  aerosol  num¬ 
ber  density  and  size  distribution  increases  with  decreasing  wavelength. 

An  increase  in  the  scattering  cross  section  results  in  an  increase  in  the 
number  cf  collisions  required  before  the  contribution  from  higher  order 
collisions  become  negligible.  For  a  source  wavelength  of  0.3u,  the  num¬ 
ber  of  collisions  required  is  approximately  20  for  the  Elterman  atmosphere. 
As  the  source  wavelength  increased  to  0.,  7y,  the  number  of  collisions  re¬ 
quired  reduces  to  approximately  5.  For  a  wavelength  of  0.7p,  single  scat¬ 
tering  is  a  significant  part  of  the  total  scattered  intensity.  Therefore, 
the  magnitude  of  the  total  scattered  intensity  is  highly  dependent  on  the 
shape  of  the  aerosol  phase  function.  As  the  source  wavelength  decreases, 
the  importance  of  multiple  scattering  increases  and  the  importance  of  the 
single  collision  intensity  decreases.  From  an  examination  of  Figures  38 
and  39  it  can  be  seen  that  with  an  increase  in  the  importance  of  multiple 
scattering,  the  dependence  of  the  total  scattered  intensity  on  the  shape 
of  the  aerosol  phase  function  decreases. 

The  results  obtained  from  Monte  Carlo  problems  in  which  the  phase 
function  for  Rayleigh  scattering  was  used  for  aerosol  scattering  showed 
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differences  of  as  much  as  a  factor  of  two  in  the  total  attenuation  for  a 
given  ground  range  when  compared  with  the  Monte  Carlo  results  for  four 
different  aerosol  phase  functions.  The  aerosol  phase  functions  used  in 
these  problems  were  for  size  distribution  models  that  are  used  by  many 
investigators  to  represent  the  aerosol  size  distributions  usually  found 
close  to  the  ground  level.  The  scattered  intensities  obtained  from  those 
problems  which  used  one  of  the  four  aercsol  phase  functions  were  found  to 
not  differ  by  more  than  35  percent  at  a  given  ground  range.  The  results 
also  show  that  the  error  introduced  by  using  the  Rayleigh  phase  function 
to  represent  aerosol  particle  scattering  varies  with  the  source  wavelength. 
For  a  wavelength  of  0.3u  the  error  is  about  20  percent,  at  0.45u  the  error 
is  about  35  percent,  and  at  0.7p,  the  error  is  as  large  as  a  factor  of  two. 

5.4  Dependence  on  Altitude  Variation  of  Aerosol  Coefficient 

Following  the  running  of  the  LITE-I  problems  in  which  the  altitude 
variation  of  the  aerosol  scattering  coefficient  was  taken  to  be  that  given 
by  Elterman's  model  for  a  "clear  standard  atmosphere,"  additional  data  on 
the  altitude  variation  of  the  aerosol  scattering  coefficient  became  avail¬ 
able  from  an  analysis  of  the  results  of  a  searchlight  experiment  reported 
by  Elterman  (Reference  4).  Elterman's  analysis  of  the  searchlight  experi¬ 
mental  data  showed  a  much  larger  aerosol  scattering  cross  section  at  alti¬ 
tudes  above  ^4  km  than  that  given  by  the  "clear  standard  atmosphere"  model 
(see  Figure  6).  The  use  of  the  aerosol  scattering  coefficient  obtained 
from  the  searchlight  experiment  in  LITE-I  calculations  would  be  expected 
to  have  a  significant  effect  on  the  magnitude  of  both  the  scattered  and 
direct  intensities  for  point  sources  and  receivers  located  above  5  km. 
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Changes  In  the  magnitude  of  the  aerosol  scattering  coefficient  at  alti¬ 
tudes  above  5  km  should  not  effect  the  scattered  intensities  when  the  source 
and  receiver  are  both  located  below  1  km  as  was  the  case  for  the  problems 
run  on  LITE-I  that  are  tabulated  in  Table  IX.  To  check  this  assumption 
the  problems  which  were  run  for  0.3,  0.45,  0.5,  0.65,  and  0.7p  point  iso¬ 
tropic  sources  in  the  Llterman  "clear  standard  atmosphere"  using  the  Haze 
C  phase  function  for  aerosol  scattering  were  rerun  using  the  altitude  dis¬ 
tribution  of  the  aerosol  scattering  coefficient  as  obtained  from  the  search¬ 
light  experiment.  The  results  obtained  from  these  LITE  problems  were  found 
to  not  differ  by  more  than  +  10  percent  from  the  scattered  intensities 
given  in  Table  IX  for  the  Haze  C  phase  function.  These  results  indicated 
that  when  the  source  is  located  at  an  altitude  of  0.8  km,  differences  in 
the  altitude  variation  of  the  aerosol  scattering  coefficient  above  £4  km  have 
no  significant  effect  on  the  scattered  intensities  at  receiver  positions 
near  the  ground  surface. 

5.5  Dependence  on  Cloud  Altitude 

To  study  the  effect  of  cloud  altitude  of  the  scattered  light  intensity 
as  a  function  of  the  source-receiver  range,  several  LITE-I  problems  were 
run  for  a  0.55  micron  point  isotropic  source  located  at  an  altitude  of  807.7 
meters.  Receivers  were  placed  at  horizontal  ranges  of  4.023,  12.07,  20.12, 
32.19,  and  48.28  km  and  at  an  altitude  of  6.096  meters.  The  ground  surface 
was  treated  as  a  Lambert  reflector.  The  bottom  of  the  cloud  was  treated  as 
a  reflection  surface  with  a  number  albedo  given  by  the  expression 

a(6  )  ■  0,639  -  0.154cos9 
o  o 
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The  polar  angular  distribution  of  the  number  albedo  for  the  cloud  was  assumed 
to  be  a  cosine  distribution.  Although  the  number  albedo  used  in  the  calcula¬ 
tions  for  the  cloud  came  from  the  Monte  Carlo  calculations  for  0,45  micron 
light,  it  was  felt  that  the  same  value  would  be  applicable  for  0.55  micron 
light  since  these  were  only  small  differences  in  the  cloud  extinction  cross 

section  for  these  wavelengths  and  the  shape  of  the  distribution  function. 

0 

P(8,A)  -  2v  j  F(9',A)sin0'de’ 

0 

for  A  -  0  45  was  a  reasonably  accurace  approximation  of  the  distribution 
function  for  A  *  0.55,  In  the  above  integral  the  function  F(9,A)  is  the 
normalized  phase  function  for  volume  scattering. 

The  magnitude  of  the  aerosol  cross  section  as  given  by  Elterman  for  a 
"clear  standard  atmosphere"  was  modified  so  as  to  give  a  value  of  the  extinc¬ 
tion  cross  secticn  of  0177  km  ^  for  extinction  between  altitudes  of  0  and 
800  meters.  The  phase  function  used  for  aerosol  scattering  within  the  atmos¬ 
phere  was  taken  from  measurements  made  at  Camp  Haven,  Wisconsin  on  31  March 
1946. 

The  results  of  the  calculations  giving  the  scattered  intensities  for 
cloud  bo* lOdi  altitudes  of  1,  2,5,  and  5  km  are  presented  in  Figures  40  and 
41  for  ground  albedos  of  0,0  and  0.9.  It  is  noted  that  for  both  values 
of  the  ground  albedo  the  scattered  intensity  at  a  ground  range  of  4  km  de¬ 
creases  with  an  increase  in  the  altitude  of  the  bottom  of  the  cloud.  For 
horizontal  ranges  of  12  km  and  greater,  the  scattered  intensity  increases 
with  an  increase  in  the  altitude  of  the  bottom  of  the  cloud. 
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CLOUD  BOTTOM  HEIGHT  (km) 

Fig.  40.  Dependence  of  the  Scattered  Intensities  from  a  Point  Isotropic  0.55/<  Source  on  Cloud  Bottom 
Height:  Ground  Albedo  =  0.0,  Source  Altitude  =  807.7  m,  Receiver  Altitude  6.096  m 


SCATTERED  INTENSITY  {photons  m_J/source  photon) 
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Fig.  41 .  Dependence  of  the  Scattered  Intensities  from  a  Point  Isotropic  0.55//  Source  on  Cloud  Bottom 
Height:  Ground  Albedo  0.9,  Sou  ce  Altitude  77  807.7  m,  Receiver  Altitude  6.096  m 
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The  effect  of  changes  in  the  magnitude  of  the  ground  albedo  on  the 
scattered  intensities  varies  with  the  altitude  of  the  bottom  of  the  cloud. 
For  a  cloud  bottom  altitude  of  5  km,  the  scattered  intensities  at  all  re¬ 
ceivers  for  a  ground  albedo  of  0.9  is  a  factor  of  about  2.25  greater  than 
that  computed  for  a  ground  albedo  of  0.0.  As  the  cloud  bottom  altitude  is 
decreased  to  an  altitude  of  1  km,  the  ratio  of  the  scattered  intensities 
for  a  ground  albedo  of  0.9  to  that  for  a  ground  albedo  of  0.0  increases 
to  a  value  of  about  3.7. 

The  following  values  were  computed  for  the  direct  intensities  at  the 
indicated  horizontal  ranges: 

Horizontal  Range  Direct  Intensity 

(km)  (photons  m*^source  photon) 

4.023 

12.07 

20.12 

32.19 

48.28 

To  check  the  accuracy  of  the  cloud  reflection  calculations,  a  LITE-I 
problem  was  run  for  a  1  km  cloud  bottom  altitude  in  which  the  cloud  was 
treated  as  a  scattering  medium  instead  of  a  reflection  surface.  The  re¬ 
sults  of  the  problem  were  found  to  be  in  good  agreement  with  that  obtained 
from  the  LITE-I  problem  in  which  the  cloud  was  treated  as  a  reflection 
surface.  Treating  the  cloud  as  a  reflecting  surface  in  the  LiTE-I  calcu¬ 
lations  required  25  percent  less  machine  time  than  did  the  problem  in 
which  the  cloud  was  treated  as  a  scattering  medium.  Thus,  it  is  seen  that 
a  considerable  savings  in  machine  time  would  result  by  treating  the  cloud 


3.568x10 
2. 369x10" 10 
4.921x10" 11 
8. 39 8x10" 12 
1. 236xl0~12 


as  a  reflection  surface. 
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